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NOTE 


Several parts of the thesis have been published 
or have been accepted for publication, 

(a) Chapter 3 (.Section 3,3) has been accepted 
for publication in Journal of Physics F 
(avv’aiting for publication) under the title, 

"An Augmented Space Formalism for Electronic 
Density of States in Cu^^^Ni^ Alloy..", 

(b) Chapter 4 (Section 4,3,5) has accepted 
for publication in Journal Physics C; 

(1985), 4677-4789 under the title, "Electrical 
Conductivity in Random Alloys", 

a) In Chapter 3 v/e have reviewed the old CPA 
theories (Kirkpatrick, Velicky and Ehrenreich- 
1970; Stocks Williams & Faulkner - 1971) 

of density of states in random alloys to 
point out the fine structures in the- 
impurity subband in density of states within 

2CPA for the same system, 

b) In Chapter 4 we have reviewed' the CPA 

work in static conductivity by Velidky 

( 1969) and Brouers and Vedyayev, 

(1972), Niizeki (1977), Niizeki and Hoshino 
(1979) to give a consistent picture regarding 
the developmoit of the theory of conductivity 
in random alloys. 
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SYN’OPSIS 

The problems of electronic states and electronic 
transport in disordered alloys are a matter of interest 
for a long time both from theoretical and experimental 
point of view. Basically the electron or lattice wave 
propagation in a disordered solid is considered as a succ- 
ession of elem.entary scatterings on the random atomic 
scatterers which are averaged over all configurations of 
atoms. For a real system in principle we have infinitely 
large number of configurations and so direct averaging is 
almost an impossible task in reality, Taylor and Soven 
(1967) first solved the problem by transforming the system 
of random scatterers by an effective medium through self- 
consistent choice of an effective scatterer which replaces 
each of the random scatterers in the frame-work of multiple 
scattering theory. The effective medium is determined 
corresponding to the choice that the single scatterer embedded 
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therein should produce no extra scattering on the average. 

The site eneroy corresponding to this effective medium deter- 
mined for a site diagonal Hamiltonian, in the frame-work of 
single site scattering theory is known as co-herent 
potential ( CP) , 

The co-herent Potential Approximation (CPA) neglects 
the scattering from statistically coupled clusters of random 
scatterers In situation of strong scattering where off- 
diagonal disorder is strong or both the diagonal and off- 
diagonal disorder are strong the need to go beyond single 
site CPA is very crucial. The difficulty in self-consistent 
generalisation of CPA suitable for describing multiple scatter- 
ing from a cluster of random scatterers has been recognised 
for sometime. The difficulty was two-fold firstly the 
summation of appropriate scattering diagrams from various 
sizes of clusters is very cumbersome and secondly the simplified 
approximation made to make the problem tractable prior to 
configuration averaging leads to unphysical results because 
of the non-analytic character of the approximate Green Function, 
The self-consistent appi'oximation appropriate for describing 
cluster effects known as cluster Co-herent Potential Approxi- 
mation (CCPA) in the frame-work of Augmented Space Formalism (ASF 
introduced by Hookerjee ( 1973) ha.s been very successful for 
various realistic model systems. 

The Augmented Space Formalism (ASF) ensures the 
essential Herglotz Property of Green Function so that the 
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reality of eigen energies, positivity and single valuedness 
of density of states are retained. This Formalism has been 
discussed in detail in Chapter 2, ASF coupled with the 
Recursion Method introduced by Haydock, Heine and Kelly (1972) 
form the practical basis for calculating density of states. 
The results of self-consistent cluster calculation of a single 
orbital model on the diamond and.b.c.c, lattices have been 
reported earlier. Calculation of electronic density of states 
of valence bands in OCPA via ASF has been also performed 
previously for random III-V semiconducting alloys, e.g, 
Ga^In^_^As (V,K, Srivastava, Ph.D, Thesis, 1982), In this 
work we report the calculation of density of states for 
metallic systems Cu^_^Ni^ here in Chapter 3. The special 
features that emerge because of effects of cluster and off- 
diagonal disorder have been discussed in detail also in 
Chapter 3, It is argued that the impurity band has got sharp 
structures compared to the CPA which is true for all concen- 
trations. The basic inputs are band parameters appropriate 
for Cu and Ni which are taken from reliable band structure 
calculations and the Recursion coefficient generated for 3d 
f.c.c, lattice of finite size configuration. 

Regarding the electronic transport, residual 
conductivity is a matter of primary interest for disordered 
solid. Electrical conductivity for a random alloy has been 
formulated by B, Velicky (1969) in the single-site CPA based 
on the Kubo-Greenwood Formalism. Brouers and Vedyayev (1972) 



have made application of this to a s-d system for semielliptic 
bands. Generalisation of caiductivity by incorporating off- 
diagonal disorder has been done by Niizeki (1976,1977) 

but in the spirit of CPA. The aim of our present work is to 
develop an explicit formalism for calculating the conductivity 
in CCPA within ASF, Just as the Herglotz property is essential 
for any physically meaningful approximation for the one particle 
propagator for any approximation for the response function the 
vertex correction must satisfy Ward Identity so that the 
averaged macroscopic conservation laws are not violated. Bethe- 
Salpeter equation together with Ward Identity with both diagonal 
and off-diagonal self-energies form basis of calculation of' 
conductivity. The conductivity formulation in CCPA via ASF 
is discussed in detail in Chapter 4, 

The conductivity results for a model simple cubic 
lattice with strong off-diagonal disorder is shown in Chapter 4 
to be very sensitive to the randoimiess in current and the vertex 
correction. Finally v;e have also reported the conductivity 
calculation for Cu>|_^Ni^ for realistic band parameters of Cu 
and Ni taken from reliable band structure calculations. For 
Cu-|_^Ni^ vertex correction is important but current correction 
is not appreciable because the widths of Cu and Ni d-bands are 
very close. 



CHAPTER 1 


INTRODUCTION 

During the period of the late nineteen sixties and 
seventies, the problem of electronic states in disordered 
metallic alloys received attention of solid state theore- 
ticians with the application and extension of the pioneering 
work of Coherent Potential Approximation (CPA) by Seven 
(1967). Kirkpatrick et al. (1 970), Stocks et al, (1971) etc,. 
Soven first gave the single- site CPA theory with the idea 
of introducing complex site energy (or potential) which 
characterises the effective medium in a self-consistent 
manner within multiple scattering formalism. Later on 
Velicky (1969), and Brouers and Vedyaye^ (1972) carried out 
the theory of electrical conductivity in substitutionally 
disordered systems within single-site CPA starting from 
Kubo -Greenwood conductivity expression (1957, 1958),' 

The validity of single-site CPA theory has been 
recognised through application to various model systems 
as well as to realistic model Hamiltonians suitable for 
disordered metallic alloys, Kirkpatrick et al, (1970) made 
application to Ni-rich Ni. Cu and Stocks et al (1971) did 

I Ml ^ O 

the same calculation for Cu-rich Cu- «Ni alloys for 
studying electronic density of states. Their results of 
electronic density of states was in good agreement with 



the experimental density of states results by Seib and 
Spicer (1970) and Htffner et al ( 1973}*. 

Brouers and Vedyayey; (1972) extended the theory of 
conductivity by Velicky in the spirit of CPA to a model 
s~d Hamiltonian with a view to study the effect of s-d 
hybridisation on the residual resistivity. Their model 
calculations show the qualitative agreement of Mott’s 
prediction of the deviation of Nordheim's rule due to the 
strong s~d hybridisation effects. In this way, single^ 
site CPA theory was applied to study the nature of the 
electronic states and transport properties for disordered 
metallic alloys like Ag^^^Pd^, etc, for tight 

binding Hamiltcnians, 

Ehrenreich and Schwartz( 1 976) remark that the CPA, 
is a mean field theory analogous to the Weiss theory of 
magnetism, the random phase approximation ( RPA) for 
describing Coulomb interactions in many electron system, 
or the Vander Waals theory of liquid gas system. While 
these theories are rigorously valid for only a limited 
range of characteristic parameters (for example, very high 
electron densities in the case of the RPA), they often 
provide physically reasonable and, indeed, semiquantitative 
descriptions in regimes substantially outside their region 
of validity. The same is true for the CPA, 

Single-site CPA theory has the following drawbacks: 

(i) it takes exactly into account only single-site 



scattering and correlated scattering from 
cluster is not accounted for, 

(ii) the effect of potential fluctuations is auto- 
matically suppressed while doing such an 
average 

(iii) it can not properly take into account the 

off-diagonal disorder in the Hamiltonian, and 

(iv) the vertex correction (i,G,, back-scattering 

effects) vanishes as an artifact of short— ranged 
potentials in single-site CPA, 

The most obvious extension of the theory to over- 
come these limitations is to go beyond the single-site 
approximation by doing a general formalism so that the 
effects of scattering fron statistically coupled cluster 
of various sites and that of the off-diagonal disorder 
could be included. Nickel 8; KruaihansJ. (1'971) tried t a make such 
generalisation using a method based on the corrected 
cumulant expansion to find electronic density of states- • 
of an one- dimensional model of an alloy, Butler (1972 , 1973) 
also tried to calculate both the real and imaginary parts 
of the Green Function using a self-consistent cluster 
method. We will not try to describe these methods here, 
except to say that they are both' Intended to improve in 
the CPA by going beyond the single-site approximation 
through making application to one-dimensional tight binding 
models. Such a genefalisation poseci severe analytic problems 
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such as having discontinuities at various energies but even 
worse, they led to multi-valued density of states. Nickel 
and Butler (1973) speculated that since this non-analyticity 
had occurred in these calculations it might occur in many 
others. 


It was pointed out by Haydock et al, (1972),. 

Mookerjee (1973), MuLler-Hartman (1973), Mills & Ratanavararaksa 
/( 1 978) 

/.that for any real potential function the Green Function 
must have certain mathematical properties which are summed 
up by the statement that it is HergJotz, Hwever, an 
approximation to the average Green Function will not be 
Herglotz, if it is obtained by ignoring contributions to 
it without regard to the preservation of this property. 

So if one makes such a level of approximation to Green 
Functions that Herglotzicity is lost, the nature Of Green 
Function may lead to unphysical results like multivalued 
density of states or negative density of states etc. In 
order to get physical results Herglotzicity should be 
retained in any approximation to the Green Function, 

The fundamental difficulty with most theories that were 
proposed during the late 1960s and early 1970s as extensions 
of the CPA is that they lead to approximate Green Functions 
which are not Herglotz, It is now realised that the preser- 
vation of the Herglotz property must be one of the central 
conside la tions in the development of a theory v\hich includes 
scattering from clusters, short— ranged order, off-diagonal 
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disorder, or the effects of positional disorder^ 

There have been attempts to include some of the 
effects of scattering frcm clusters of atoms, atleast 
conceptually by using a technique called the molecular 
coi-herent potential approximation (MCPA) as discussed 
by Tsukada (1969), and Butler (1973). In fact, it is a 
straightforward generalisation of the CPA approach in which 
the lattice sites of the alloy are partitioned into cells 
that contain a number of sites, which is greater than one. 
Mathematically, it is obvious generalisation of the 
T-'^roatrix equation in CPA, The effective self-consistent 
molecular scatterer is found by solving this T-matrix 
equations and used to calculate a Green Function, One 
can prove that the MCPA leads to a Herglotz Green Function 
within the same philosophy as that the CPA does. The main 
objection why MCPA is not a satisfactory theory is the 
imposition of fictitious unit cell boundaries and as a matter 
of fact the translational invariance is lost although it 
is invariant under cluster translations. 

An approximation scheme for including the scattering 

from clusters that is very much easier to do calculations 

than the iMCPA have been! also proposed. This is known as 

embedded cluster method (ECM) as discussed by Inglesfield 

(1981), Baraff et al (1986) and references therein. The 

idea of ECM has been used by Gonis et al (1984) in Ag^ Pd 

1 — c c 

where they embedded the cluster in a CPA medium. The main 



objection to the ECM is that ”it is just a cobbled up ■'tbe-o-yy 
in the sens.e that it is not a self-consistent theoby like 
the CPA or even the I'l/iCPA” (Faulkner 1982), 

All the above techniques described do not provide a 
reliable solution to the theoretical problem of developing 
a self-consistent theory of the electronic states of an 
alloy that goes beyond the single-site approximation. The 
problem of going beyond single-site CPA theory has been 
proposed by Mills (1978) in his original studies (known as 
Travelling Cluster Approximation) on the problem of a 
Herglotz average Green Function within the frame work of 
diagrammatic perturbation theory. But the problem is very 
cumbersome to evaluate the contributions of higher order 
terms for various sizes of cluster for a real three-dimensional 
lattice. Mills and Ratanavararaksa (1978) reported the 
density of states results for the three-dimensional (simple 
cubic) tight binding model of an alloy having diagonal 
disorder only in tv\/o-site travelling cluster approximaticari 
(TCA). 

Mookerjee (I973a,b, 1975a, b,c) introduced a general 
formalism known as Augmented Space Formalism (ASF) to 
calculate configuration average of a general function of 
randan variables by using abstract operator formulation of 
quantum mechanics. He was able to establish the idea of 
incorporating scattering effects frcmm statistically coupled 
sites in a finite size cluster as well as the importance of 
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off-diagonal disorder in a self-consistent manner. 

The formalism came out with success in producing 
Herglotz analytic Green Functions through applications to 
various model systems (Kaplan and Gray 1976, 1977, 1978j 
Kumar et al* 1982; V,K, Srivastava, Ph.D, Thesis, 1982} 

Kaplan et al, 1980; Thakur et al, 1986), Here, we, for 
,the first time made application of ASF for the calculation : 
of density of states of a disordered metallic alloy like 
Cu-rich Cu^_^Ni^ corresponding to the realistic choice of 
tight binding parameters from reliable band structure 
calculations for its constituents and using Recursion Method 
(Haydock, Heine and Kelly, 1972, 1975) to generate the 
local Green Function for d-bands. Results of density of 
states have been reported here in two-site CPA in ASF which 
shows sharp structures in the impurity band. Note that 
the fine structure in the cluster CPA (CCPA) density of 
states is a ccmmon feature which has been reported earlier 
by Kumar et al, (1982) in semiconducting alloys. The 
cluster CPA within ASF has been also applied for the calcu- 
lation of vibrational properties of random transition metal 
alloys, e,g, Ni^_^Pt^, Ni^^^Pd^, Ni^__^Cr^ by R.P. Singh (1982), 
Mookerjee and Singh (1985), 

So far ASF has been proved to be successful in 
dealing with .either the calculations of electronic density 
of states or the vibrational density of states. It is a 
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general formalism in the sense that it can be applied to 
calculate configuration average of any physical observable 
which is a function of random variable provided that 
probability distribution corresponding to each random 
variable is known. 

Another very important physical quantity of interest 
is the static residual conductivity. The theory of electrical 
conductivity in the strong scattering regime (beyond the 
Boltzmann limit) has been carried out by Roth (1975) and 
Singh & Roth (1982) for liquid metallic alloys in the 
effective medium approximation (EMA), The importance of 
back scattering effect has been recognised by them which 
was missed by Asano and Yonezawa (1980), The old CPA 
theories (Velicky 1969; Brouers and Vedyayeb 1972) can not 
take into account the back scattering effects because of 
short-ranged potential used in the standard single-site 
CPA which leads to zero back— scattering effect on the 
average. Niizeki (1977a,b,c); Niizeki and Hoshino (1977) 
emphasized the importance of vertex correction (which 
characterizes the ■ back-scattering effects) for long range 
potentials within a bond CPA type approach. They showed 
that static conductivity is very sensitive to the randomness 
of the current and the effect of vertex correction* The 
importance of vertex correction is discussed for multiple 
bands by Chitnavis and Raghavan (1983)# Chitnavis & Leath (1983) 
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Following this we have developed a cluster CPA 

formulation within ASF by incorporating off-diagonal 

* ' 

disorder and cluster effects from the statistically coupled 
sites for the evaluation of the configuration average of 
the Kubo-Greenwood conductivity expression. We have 
derived analytic representations for vertex correction as 
a function of CCPA Green Function§and self-energies. We 
have reported here the static conductivity results as an 
energy functional for a model simple cubic lattice with 
strong off-diagonal disorder. Finally, we also carried out 
the calculation of static conductivity as an energy functional 
in noble/ transition metal alloy like Cu^__^Ni^ within two- 
site CPA, 

In both the two cases it has been seen that vertex 
correction is not zero and the effect of randomness in 
current is appreciable for strong off-diagonal disorder. 



CHAPTER 2 


FORIVIALISM 


2.1 Augmented Space Formalism 

2.1.1 The Augmented Space Methods 

The electron in a disordered solid experiences random 
potentials at the lattice positions and the system can have 
infinitely large number of configure tions , The disorder 
in a solid is characterised by a set of hamiltonian para- 
meters which are random variables and obey some probability 
distribution. So description of such systems should be 
from a statistical point of view and one hopes to describe 
the physical properties of the system through configuration 
averaging of the physical observables. 

The augmented space formalism is an attempt to 
systematically take into account the configuration fluctu- 
ations while studying configuration averaging. The formalism 
is exact but for practical calculations approximations 
preserving constraints of physical origin may be generated. 
Augmented space formalism (ASF) was first introduced by 
Mookerjee (1973a, 1973b) and later on Kaplan and Gray 
(1976, 1977), Kaplan et al (1980) provided detailed exposi- 
tion of this. The formalism has been proved to be powerful 
and efficient method for averaging quantities like one 
particle Green Function, Response Functions related to 
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two-particle Green Functions particularly when we want to go 
beyond the single site coherent potential approximation (CPA) 
and include effects of clusters, large off-diagonal disorder, 
short ranged order and so on. Previously ASF has been applied 
to average only Green Function as a single function of random 
variables. Here, our aim is to apply ASF to average physical 
quantity like static residual conductivity which involve 
products of moi-e than two functions of random variables where 
direct decoupling is impossible, 

2,1,2 Mathematical Form of ASF and its Application to 

Model Systems; 

The detailed mathematics of ASF is carried out in 
various communications (Mookerjee - 1973; 1975a^,b,c) 

Here we briefly discuss the main points related to the aug- 
mented space formalism to substitutional disordered system. 

In a substitutional system the lattice sites are 
randomly occupied by the atoms of the type A,B,C etc. So 
the site energy e^^ is a random variable which can take on 
value e^;, eg etc. The random variable will obey some 
probability distribution p-Ce-) such that p.(e. ) > 0 and 

.+00 1 i 1 1 

- 00 ~ "that the probability 

density is analogous to density of states of a solid in 
the sanse that it is positive definite and while on integration 
should yield the number of electrons which can be normalised. 
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Mookerjee (1973-) suggested that one can introduce 
a hypothetical space corresponding to a random variable 
such that a basis set { I 1 and an operator in that 
v/ould characterise the random distribution of e^ * Note 
that the quantities I ^ respectively 

the analog of eigenvectors, hamiltonian and density of 
states in a Hilbert space. Later on Mookerjee (1975) descri- 
bed 1 f^> as the disorder field in analogy with a fermi.on/ 
boson field and was able to give diagrammatic representation 
of it in a similar way as one does in Feynman diagrams for 
fermions or bosons* But here the problem is different 
because usually knowing the hamiltonian we want to find out 
density of states but here for a given probability distri- 
bution we have to find out an operator corresponding to 
p(ej_). Idea is the following, if 1 > is a specially 

chosen member of the orthonormal basis { lf^>} in then, 

Pi( ef) = - 1 Im <f ( el - M)"'' f f j > ( 2J ) 

j_ 

where e=e^ + i-»o'. 

The problem can be handled in the analogous way to 
what one does in the continued fraction expansion of Green 
Function (Recursion Method; Solid State Vol,35 Haydock) in 
Recursion Method which v^q would discuss in the next section* 
One can try to find a convergent continued fraction of the 
probability density p(e^) of the following form 



13 



The operator M is tridiagonal matrix with a^ , a 2 , 
etc, along the diagonal and b;|,b 2 ,,, etc,..,along the off- 


diagonal positions in some basis 


It should be noted 


that for some distributions like Lorentzian one can not get 


convergent continued fraction expansion. So the whole thing 
fails theroi Let us choose our system to be a random alloy 
^c®1-c diagonal disorder (e^) only. In the tight binding 
Anderson hamiltonian 


H =5e^p. + V , (2,3a) 

the diagonal term e^^ form a set of random variables descri- 
bed by* e^ = Ce"^ + e"^ ( 1 - 5 (2.3b) 

where =* 0» 1 foi* binary alloy. 

If one neglects the short range order, probability density 
is given by* 

p(Nj^) = c 6(N. - 1) + ( 1 - c) 5(N^) 

= _ 1 In J ) 

, N, + iti- 1 N, + i 

T) - 0 ■ ^ 
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where = c? 82 = (1-c) and b‘^ = c( 1-c) , 
So M has (2x2) matrix representation 


c c( 1-c) 

'^c(1-c) 1-c 


(2.5) 


with eigenvalues 0 and1 . Its eigenvectors Ivi> and fv^> 
describe occupancy states corresponding to the eigenvalues 
0 and 1 , Now we have to define general state vectors as 
for example f f^; r^> = If^>S:<lrj^> in Augmented product 
space W = XTtqjj^ and the average becomes, 


+ 00 /;j\ , 

F = < f^? r. I[ / dm (h + iO ) >x 

-00 

E( h) dh 3 [ f^, r^ > 
where p(e) = - (h + io"^) , 


( 2 . 6 ) 


If wo want to generalise the augmented space for 
a set . { e - } then ^ The basis in the product 

space is 

I f > = 5lf]^> {If?> .... 

and if g^^ ^ ( z) = ( zl - M^)"^ then 

P({e^}) = - ^ Im <fQ I G( z) I f^> (2.7) 

where ! f ^ > = ! fj, > I f ^ > 0 I .... 

G(z) = * Hz) « HU) 


and 



Let us consider the probability distribution for ,a 


single variable e 

_ -hoo 

f = f( e) p( e) de 

= S f( e) - Im ( z) de 

— 00 

z -* e -1- iO 

where, p( e) = - Im • (2,8a) 

Let us consider f( z) a function of a complex variable z and 
has no singularities on the real axis in the neighbourhood 
of the branch cut of the function z) • 

f = 

or 7 = 

where g ( z) 
operator of M, 

Thus f 

= (2,9) 

where f(M) is the same function of M as f is of e. So we 
see that in general the average can be represented as a function 


^ « f(z) dz 


+ 00 

'I f(z) dz <f 1 ( z-h)'"'' dp(h)lf> 


2Tti 


-oo 


o 

(2.8b) 


"HCX) t /. , \ 

= 7 snd p(h) is spectral projection 


--00 


z h 


H'OO 

■<f I ^ f(h) dp(h) I f > 
-00 



I D 


Of the constructed operator M which involves the probability 
distribution. 


2.1,3 Hamiltonian in ASF for Both Diagonal and Off-Diagonal 
Disorder; 

Our starting hamiltonian is of ti(^'it-binding type 


» = ^ "in hn + 2 

1 n n,m 


V T 

in,jm in,jni 


where e^^j^ 


A 

^in "i 


C N, + e®^ 


(2. 10a) 
(2,1 Ob) 


i,j being site-indices and m,n band indices. 


V. . = N.N. + V®® (1-N. )(1*"N.) + 

in,jm nm i j nm '' i'^'^ j-' nm 


iN^d-Nj) + N. (1-N^)1 


(2,10c) 


Here e- represents diagonal disorder and V. 

in f j ^ 

represents off-diagonal disorder arising from site energy and 
hopping integrals for the constituents A,B respectively. 


p(Ni) 


H 


= c 6(Nj,-l) + (1-c) 6 (N^) jN. = 0,1 (2.11) 

_ _ ■ ( 1 ) 

I u Ij \r N,N. 

“ ^B i n . n. ‘ i^in ^ i^tj n,ra ^nm i j 




( 2 . 12 ) 


T 

xrifjm 


S 
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where e„ = ( - e^), + vf ^ - 2V^® and 


nm nm 


v(2) = v-'® - V®® 

nm ''nm ''nm 


(2.13) 


As before M = 


C ''‘cd-c) 

c( 1-c) 1 -c 


in the basis IVq> ,lv^>, 

(2,14) 

The augmented space hamiltonian is 


H = Hd fx'' I + p e 
® ' n,i 


Pin @ ,x) 


.^.^nm^ ■^in In ^ <x-', (5c) + 

nm in^jm 


2 . 5 . T. . 5 f (m(^^ f 3 fi dC' 

i.m nm xn.nm '4’ L^-’ 'X - 1 - 


n,m ’nm ij^j ‘in,jm 

I 1 being identity operators in all subspaces jlr (2,15) 

except to those superscripted. 

Until now we have not made any approximation. But the evalu- 
ation of Green Function in augmented space is very difficult 
because the large number of self-avoiding path which starts 
from some point and comesback to the same point. This problem 
is solved by Graphical Methods in conjunction v\/ith suitable 
forms of Green Function by the continued fraction represen- 
tation of the Green Function through Recursion Method, 


The Augmented space formalism is not just a mathe- 
matical way of averaging but one can analyse the cluster 
effects coming from statistical coupling of various components 
of the system. Before going to discuss its physical relevance 
to cluster co-herent Potential Approximation (CCPA) one 
should briefly discuss the Recursion and Graphical Methods 
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because we are going to use the to calculate density of 
states and conductivity of the system. 


2 ,2' Graphical Meth od 
2,2,1 Introduction; 

In this section, by a graph we shall mean a set of 

'vertices’ connected by ’links'. The Graphical Method was 

first introduced by Haydock (1972) for calculation of resolvent. 

We would be using it to evaluate Green Function matrix element 

in augmented space and one sees the approximation made therein 

is related to cluster co-herent Potential Approximation, A 

graph is associated with every basis in { Ii> 1, To every 

element ! i> of the basis set {(i>) is associated a vertex 

and to each distinct pair of elements ( I i >, I j >) a link 

or bond 1. . which may be directional (i,e, 1 . . ^ Suppose 

we want to invert the operator M then the contribution of each 

1 

vertix v^^, is defined as k(v^) = 


,, and that of the link 1. . 


tv. 1,. ^•v,. ,,,. Vj.j 3, Contribution of a path is defined 


as a series of ordered vertices and links, i,e,, 
as kCQ„) 

^vertices 

‘Now forG = (zl~ H)” ' the matrix element 


% 


k(l,,) . 


links ij 
^-1 


G... = < i 1 (z I - H)”^! j > 




00 ^ 

Nio 


( 2 . 16 ) 


where are all self-avoiding paths from v^ to v^ , By a 
self-avoiding path we mean it does hot trace any vertex more 



than once Inbetween two points for a given path. 


2,2,2 CPA by Graphical Methods in Augmented Space; 

In Augmented space the basis set { I nf?> } contains 
two informations spatial extent of the lattice by n and 
{f^the disorder field information. Now for a given probability 
density p(e) of site energies of lattice points and conside 3 >- 
ing only diagonal disorder in the hamiltonian one can write. 


-nf,mf« 

P( e^) 

I f > 


«nm . 


ff 


^ Im <f [ (2 I - M^r'' I f > 


! fl > (x; I f^>(x) 


f being the ground state. 


(2,17) 


(2,18) 


An electron at a site labelled by n and field state__^ 
If£> can either make spatial hop to one of the near neigh- 
bours of the nth site with matrix element V keeping the field 
state same or it can remain in the spatial site while the 
field at the site n changes according as IV?" , the field at 
other spatial sites remaining the same, Suppose we want to 
find 


Gnm( z) = < nf I ( 2 I - H)"^ I mf > 

The graph is shown in Figure (2 ,1a) for G^^ 


(2.19) 
and for 
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For an ordered system we have 







( 2 . 20 ) 


where Rq( z) 


^ S p(o)^ ^,3 are neighbours of zero and 

Rq( z) is the contribution of all nonintersecting 
avoiding paths from vertex 0 and back in the spa 
In the augmented space the configuration average 
Function G( z) can be written as. 



( 2.21 ) 


R( z) is the contribution of all self-avoiding path from the 

f \ "^*1 T'^lv m the 

vertex of and back that are either (a; 

oelf-avoiding paths 

spatial part of the augmented space or \.dj 
in the augmented space which Include field hops b 

foim closed loops. i z) is the contribution of all closed 

1 -vn +hp full augmented 
self-avoiding loops fron Of and back m 

i.. 1 via field hops. Now 

space connecting different spatial sir 

in single site CPA we delink all T( z) contributions. 


Under this approximation, we have. 


G°Hz) = 


[2 - Ro( z - *<t 
Poo - "1 


_ tv^'ajTz)3 


( 2 . 22 ) 


which is clear from the graph shown in Figure (2.1f) 
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We have also, 

z) = — ; ^ 

z - Rg (z - W^G, ) - Si+i(2) 

(2.23) 

Now, in the conventional CPA, we have our self-energy 
2 ( z) defined by 

G (z) = — — — - 

[ z -r ( z) - Rq( z-2) J 

= Poo*^*^^ (2.24) 

Equations (2,24) and (2,22) show the equivalence,. Equivalence 
of the two has been discussed elaborately by Mookerjee 1974, 

2,2.3 Cluster Co-herent Potential Approximation (CCPA) in 
the Augmented Space by Graphical Methods: 

We have derived earlier the representation of the 

operator in disorder fieli space that characterises 

the probability distribution p(ej^) in a basis I v^ > , I v^ > 

in equation (2,5), 

The corresponding operator representation of the 
hamiltonian in augmented space which couples the represen- 
tation of the hamiltonian in the Hilbert space to the disorder 
field is also shown in equation (2,15), 

Wo start from the configuration averaged Green 
Function, 



(l2f) (lit) 



(0,fo) 


(a ) 



© = Ob 

^ 7 Z7^ = oj 

Represents an octagonat 
decoration as In fig . 2.2 a,b 

(b) 


Fig. 2.1 (a) CPA Graph in ASF (b) Rehormalisation of 

site end bonds of the octagon ( 2 CPA graph } 
in ASF due to the infinite medium 
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'^in, jm^ z) = <r^n'r°I(zl - H)“^ ^ (2,25) 

0 ^ i 

with IT - 1 v^> , i,j being site indices and n^m being 

band indices. Note that upto this stage we have not made 
any approximation, VVe have discussed in the previous section 
that Green Function matrix element can be evaluated by graphi- 
cal method and there is' infact close connection between cluster 
CPA and the nature of the self-avoiding path that traces the 
solid from one point to the other in augmented space. In 
single-site CPA we have delinked all closed self-avoiding 
paths involving both spatial and field hops in the augmented 
space. 

In a (n) cluster CPA only those self-avoiding non- 
intersecting paths will be taken into consideration which 
involve disorder field hops over the n-sites of the cluster. 

For a CCPA with ( > r2»,,.rj^ ) within the cluster all 

polygonal paths that do not completely lie in the subspace 
dH and involve those vertices and none other must also be 
exactly accounts for. These correspond to the multiple 
scattering within the cluster. Larger polygonal paths 
which involve sites within and without the cluster but 
which do not completely lie in the subspace must be 
delinked and infact this is the approximation one makes 
to evaluate the configuration average of Green Function, 
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2.2,4 The Self-consistent Effective Medium in Augmented 

Space; 

V)Je have already discussed how in graphical theory 
one considers the CCPA corresponding to a given size cluster. 

So in graph theory one would be able to determine the self- 
avoiding polygonal path corresponding to the size of the 
cluster concerned. The polygonal self-avoiding path which 
takes into account the multiple scattering within the 
cluster is infact immersed in the rest of the lattice and 
so the vertices and bonds (links) should be renoimalised 
accordingly. Let us illustrate this by doing 2CPA for a 
random binary alloy. Corresponding to the 2-site cluster 
we have an octagonal decoration of self-avoiding polygonal 
path in the augmented space graph as shown in Figure ( 2,1 a, 2. 2a, 
The eight vertices of the octagon corresponds to the eight 
different configurations of the two— site cluster (2x2 ), 

The octagon belonging to a bond has no connection to 
another bond through the augmented space link. What it 
means is that the correlated scattering from three or more 
sites belonging to different bonds is ignored. 

Finally we renormalize the vertices and bonds in 
the following way. The renormalization is a two-step 
process. 

Step 1 ° 

The decorating octagon is not an isolated octagon 
involving only the two spatial sites (O and 1), If we look 
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at the vertex A ( 1 f ^ ) and the bond AE ( 1 f ^ to Of^ ) belong- 
ing to the octagon is only one of the z bonds emanating 
from A to the z nearest neighbours of 1 , The other remain- 
ing (z-1) bonds (from 1 f .j to 2f ^ , 3f^,...etc.) also hang 
on the site A. This is also true for other vertices (G, 

H, etc.) as well. So the whole octagon has to be renor- 
malized due to the effect of the rest of the lattice in 
the augmented space representation. Once this renorma- 
lization is. carried out let the whole octagon be now iso- 
lated. Mathematically the renormalised vertices z) 
and bonds o^(z) are found as follows. 

Suppose the whole lattice is divided into two 
subspacesj (i) An unrenormalised bond AE which we call 
system ( 1 ) It has a hamiltonian 

= V2 (T^g Tg^p (2.26) 


and (ii) A lattice L which is the original lattice minus 
the link AE and in which all bonds and sites are renorma- 
lised by S . We call this to be the system (2), This has 
a hamiltonian 



So(2) 


i^A,E 


Pi 


+ 2^ (2) 




"ii 


(2,27) 


Also there exists a linking hamiltonian, 

p|(int) ^ S^(z) 2 (T^^+ T^^ + + Tg^) (2,28) 

Now Green’s ooerator corresponding to the subspace 1, 
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;( 1 ) - 


_p 


X 

G 


P. GP = [z I - _ H^nt ^(2) ^int] 

-I 


Inverse of the operator X in subspace Y 
Green's operator in. the total space ,and 

r ( 2 ) _ / T J 2)^"^2 

G = ( z I - witn 


-P. 


(2.29) 


(2.30) 

The above equation immediately tells us that the 
effect of the rest of the medium on the bond AE is to 
modify the hamiltonian by 



given by 


H 


( 1 ) _ 


H 


( 1 ) 


+ o where the self-energy ^ is 


AA 


EE 


= ^o(z) 

S { ( z) e jj,. ^km' (2,31a) 


^AE “ "^EA 

= (z) 

^ keNy^^meNg *^km (2,31b) 

= Nearest neighbour of A 
Ng = Nearest neighbour of E, 
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N'ow as stated earlier that this isolated renoimalised 
octagon has to be also renormalised because each of the 
vertices A,E has got connections with other vertices through 
either a field hop or spatial hop. This final renormali- 
sation would be carried out in the following way, Vife divide 
the octagon into two subsystems a bond B(Of, If) and the 
rest of the octagoi M, So that the system has the following 
hamiltonian representations, ( see Figure 2,1b, 2-2a,b for 
renormalisation in graphical methods), 

S “ i,je0f,1f Sij' - " ij?^0f,1f -ii » 


r>> 


2 2 (H 

iEB jeM '^"ij 


(Hij Hji) 


(2..r 

(2*3c 


Now, we have, 


G = Pb^Pb 

-p 

r ^ ~B «i^BM 1 ^ 

= [zI-H-ri GB J 


pM r 

G = Lz I 


h^3 


(2,3<^ 

(2,34b) 


Now we have the final renormalised site 2^(z) and bond 2^ (z) 
givaa by 


-M uBM 
°Jk “ko > 

(2,3? 

„BM pM „BM 
^oj ^jk ^k1 

(2.3c 
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Equations( 2.35) and ( 2.36) are the set of self consistent 

form' of equations which would provide self-consistent 

effective medium = ■2^)* 'the forthcoming 

^ o 

chapter, we would make the application of 2CPA to 
random substitutional alloy and discuss the physical effects 
of clustering coming from the statistical coupling of the 
two components of the system, 

2.3 Recursive Solution of the Schrodinger* s Equation 
2.3.1 The Recursion Methods 

The basic philosophy behind band structure like 
studies of the electronic structure in a solid is to exploit 
long range periodicity of the atomic potentials and to 
express the electronic properties of the solid as a coherent 
superposition of the electronic properties of all the atoms. 
But this is true when the interaction between electrons 
and atoms are not too strong. When the electrons interact 
strongly with the atoms this picture breaks down and prope3>- 
ties no longer depend on long range periodicity but rather 
on only the first shells of neighbours of each atom. The 
d electrons in transition metals are prime examples of this 
regime. The physics is better understood by means of a 
solution that explicitly accounts for the role of local 

environment. Underlying physics in recursion (Solid State Ph . 

1^* 

35 (1980) -• Haydock) is that the local orbital itself 
has the greatest effect and that successively more distant 
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rbital, have lesser effects* 


orbi 

The solution must therefore define a hi 

environments so that their relative influonct: 

PhvsicallYj 

displayed in the local density of st- x- • 

. -ir^+pnqitv of each eigen- 

the local density of states is tne in 

1 - j Thp local density of 
state on a particular atom or band. 

. , is determined by the sequence 

states for a basis orbixal 

, ?u 1 where 

of environments described by orthogona n 

{ 1 is a linear combination of basis orbita o 

is . The !up must satisfy the condition that H u„ = 

“n ^ “"n+l “h+1 '’n Vl ’ 

^ n ^ Cir6 

H is the hamiltonian of the model an jq 

on the Shell of atoms n hops from the atom accommodating ^ * 
The parameters a, and b^ describe the coupling of each environ^ 

ment to itself and its neighbours. The aim of g , 9 

to such basis set is that one can transform any 3 

to pseudo- one- dimensional chain model where 

and b the coupling cons tan s 
characterised by the parameters ^n "" 

of each environment to itself and its neighbours. 

XX 0 can be written 

example the local density of states o 

X e ihe chain model as the continued 

in terras of the parameters of the chai 

fraction 


h„(E) = 


, 2 // p o /( E-arj^bo'^ • * • ^ ^ 

Cim l/CE-a^ " 1 ^ 2 . 


(2,38) 
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where 

nn 

Let us try 


^ ^nt 17 n • 

to get the continued fraction for local density of 


states c^(,E) starting from the Green Function, One can 
immediately check that in the above basis set {ru^>l 


hamiltonian is a tridiagonal matrix 




^0 


0 0 

’ 






b^ 0 



H = 


0 

^2 

^2 ‘’a 





0 

0 

b^ a^ 

; 


1 

i 


• 

• 

3 3 

•• • 



! 

! 

V 

{ 


« 

• 

• 

« 

• # 

• V 



1 

Go(E) 

= 

O 

V 

r (21 

- H)"'' 

O 

V . 

I 

1 



i 

(E-a„ 

) 

-^1 

0 0 

j 

i 

i 


i 

1 

f 

-b^ 

( 

E-a^) 

~^2 

i 


1 

0 


-bg 

(£- 82 ) -^^3 



i 

i 

0 


0 


• • • • j 


i 

i 

i 

« 

0 


• 

0 

♦ • 

• • 



! 

i 

f 

m 


0 

• • 



Let us define the determinant of the matrix with the first n 
rows and n columns deleted as Dj^(E), 
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Gq( E) 


D^(E) 




(E-a^)D^(3) - D^CE) 

1 


(E-a^) - b 


2 ^2^ ^ 

1 iT^^xr 


1 

E - G,(E) 

If we continue the process n steps we have, G^, G 2 G^, 


Go(E) 


.-I . _ 

, T - ; - ^ ■ 

E — h^ XE— 3 j—b^/ • * •XE'*3j^__^"'fc*y^G^( E) ) 


(2.39) 


where, 


G„(E) 



1 

‘’n+1 Vl '('e) 


We see that G^CE) has continued fraction represen- 
tation which in general may be very large in extent. It 
has been verified that for many systems the continued 
fraction converges. When the convergence is slow or poor 
it can be terminated by some termination schemes. So 
after doing reasonable level of recursion, for example 
10 level is sometimes sufficient one uses termination scheme 
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as developed by C.M.M. Nex (1978, 1984), Haydock and Nex 
(1984) . 


2,3,2 Theory of the Continued Fraction Terminator; 

Underlying philosophy for termination of the continued 
fraction for Green Function is that the approximate resolvent 
R(2) (which replaces G) should be such that the corresponding 
hamiltonian H should have the similar energy spectrum. By 
similar energy spectrum we mean it should correspond to as 
many as possible of the same singularities as the true 
hamiltonian H should have, i.e, band edges, Vanhove singula- 
rities etc. The continued fraction form of the Green Function 
for a uniform chain parameters, i.e, a^ = a, b^ = b, 


Go(E) = 


■ — — 2 — 

E _ a - b G 


or 


b^-G^(E) - G°(E) (E-a) + 1 


Therefore, 


Go(E) 


(E-a) (E-a)^ - 4b^ 


(2.40) 


Here one sees that we need not have to use any termination 
schane because the Greai Function is exact. Now let us 
consider a general system having a single band hamiltonian. 

In this case square root terminator is appropriate as shown 
by Haydock and Nex (1984), Hers, we illustrate its application 
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msth© iTH 'ticc. lly ^ Ia/^ h^vG th 0 Gi^ggh Functioni 

‘ i 


G(E) 



- b: 


n-2 


^ “ ^n-2 ~ ^n-1 

where t(E) = Reqi^ired tenriinator. 


(2.41) 


Since an eigenstate of the chain is seme linear combination 
of the states u^, j we write the Schrodinger 

equation for the chain in the following form 


H 


00 

S 

n=0 


00 


n 


= E S 


n=0 


PnS 


( 2.42a ) 


In matrix form 


(E- 


-D, 


0 


(E-a^) 


-b, 


2 


0 

•br 


0 .... 




Pi 


n 


= 0 


(i2.42b) 
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The recursion co-efficients can be found in terms 
of orthogonal polynomials. We use those of first kind 
p^(E) and the second kind q^(E)., Now making a proper 
choice such that 

p_l = 0 = q.,, Po 

o 

and 

(E-aj)p-(E) - Pi_i(E) , 

i = 0,1, 2 .... (n-2 ) (2,.43b) 

and for i = 1 , 2 ,,.,,,.. (n-2) 

q^(E) = (E-a^)qj_^^( E) - b? qi_2^E) (2.44) 

One writes 
G(E) 

For a single band density of states 

t(E) =rE - (a+p)/2 - '^rE^TCT)]/2b^ (2,46) 


2 




n-1 


(E) - b^_^t(E) p^_2 (E) 


(2.45) 


we havep^._j_^ ( E) = 


(2,43a) 


a , -p are band— edges and 


For many band case trivial extension of this can be 
made as discussed by C,M,M. Nex(l984), The various possible 
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termination schemes has been discussed by Haydock and Nex 
( 1984 ) and .. JA.M. Nex ( 1984 ), Now-a-days Recursion Method 
IS used to find local density of states for various ordered 
and disordered systems. The newly developed termination 
schemes heve macie the method very efficient. 



CHAPTER 3 


■ PE .ELSCTRONIC DEInISITY OF STATES FOR METALLIC ALLOYS ’ 
LIEANSITION/NOBLE 1?ffALT~IN T lAL APPFiD- 

.^IJAIIQiiTCPAj AND CLUSTER CO-HERENT POTENTIAL APPRO- 
XIMTION TcCPAf 


3,1 I.ntroducti_on 

Soven (1967) first crystallised the idea of the 
co-herent Potential Approximation (CPA) for describing 
electronic density of states in disordered systems. Although 
the idea was not new. Self-consistent effective medium 
theories had been discussed in one form or another in different 
contexts? propagation of waves through disordered layers, 
dielectric response of insulators injected with conducting 
globules and so on. We should like to emphasize that the 
ideas behind the CPA and its cluster generalisations are a 
lot more general than the specific application here may 
suggest, 

CPA has been applied to various disordered alloys and 
it is reasonably in agreement with experiments. The usual 
single-site coherent poten-tial approximation arises from 
the self-consistent solution of the approximated multiple 
scattering equations. The approximation involves replacement 
of each of (N-1) scatterers out of N scatterers by an 
effective potential such that there is no extra scattering 
on the average from the Nth center representing the exact 
potential, Soven ( 1967) and Velicky et al (1968) made 
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application of CPA to various model systems. Our intention 
would be to generalise the CPA to include contributions fran 
exact clusters self-con sistently. The aim would be to apply 
the formalism to study the electronic density of states for 
disordered metallic binary noble metal/ transit! on metal 
alloys. Calculations of- density of states for Kfi rich 
Cu-Ni alloys within the CPA frame work has been discussed by 
Kirkpatrick, Velicky and Ehrenreich (1970), Later on calcu- 
lations on Cu rich Cu-Ni alloys have been shovm by G.M, Stocks, 
R,W. Williams and J.S, Faulkner (1971) through trivial 
extension of the model by Kirkpatrick et al, (1970), 


3,2.1 The CPAj A Recapitulation: 

Let us recapitulate single-site CPA theory and its 
application to noble/transition metal alloys. We woul^^, |je 
discussing the effective medium approach to derive the CPA 
equation. The well known tight binding hamiltonian for a 
single-band has the following form, 


H 


1 11 





where, for binary alloy for example, 


e 


B 


(3,1a) 






site energies of A, B 

= Hopping integral corresponding to a hop from 
ith site to jth site 

= Transfer operator; = projection operator 
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where 


‘ ) 




site indices. 


o c ^ 

More explicitly, (p /2m + V^)Ii>= I i > 


(3.1b) 


V- = crystal potential corresponding to the ith site 


and 


V. . 

ij 


-<iln^i V^U> 1 


(3.1c ) 


Suppose the hamiltonian written in equation (3,1) has only 
diagonal disorder, i,e, is only random. Then underlying 
philosophy in single-site effective medium theory is that 
if any of the effective potential^ ( z) (Co-herent Potential) 
is replaced by an exact potential there would be no extra 
scattering on the average. Mathematically this means the 
average of the single-site scattering operator 

\ (z) =(e. -S(z))/Jl _ C (3.2) 

vanishes, i.e., 

<'^[2(z)]>= 0, - (3>3) 

whe re , 

S ( z) = self-energy or Co-heront Potential 

r. , . + 

z = E + 10 

and 

F( z) = Green Function for the alloy. 

More explicitly, F( z) = F^ (z-2(z)) (3.4) 

F^ represents the Green Function corresponding to the 
crystalline system. Equation (3,2) together with equation 
(-3,4) constitute ttie self-consistent equation in single site 
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Co— herent Potential ApproxinE tion which has the following 
form? 


c {/a 


{1-c ; (e -S (z)) 

= 0 

1 *- (^B 

(3.5) 


Equation (3.5) can be written in more convenient form as 
follows; 

c(1-c) F .( 2 ) 

2 (z) = e + - — 4_ ; — 

1 _{(i_2c)6 + (£-2(2) 


)J F(z) 


(3.6) 


where, 




We have made the assumption that the hamiltonian 
involve single-band parameters £p, e b> ^AB^ • 

But this assumption is far from reality. A lattice site in 
a solid can have in general many orbitals or bands. For 
example in a system like Cu^_^Ni^, c being concentration, 
Cu— site or the Ni site have s and d bands both of which 
play significant role in its various physical properties 
magnetism, resistivity, specific heat etc. It has been 
pointed out by Hodges et al (1966) that the s-d hybridisa- 
tion plays important role in transition and noble metals. 

So one should start with the density of states of the 
pure system where hybridisation is taken care of to do any 
realistic level calculation for such random systems. 
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Our aim wou3.d be to consider the general hamiltonian 
in v^^hich hybiridisation effect is taken into account. The 
general hau'iltonian of s-d noble/ transition metal system 
has the follovd.ng form; 

H = lns.>e=<nsl + ^ ''nm ' 

- 1 - ^ f nd>E^ (n) < nd I 

+ I Y„(x) [ lns>< nd I 

+ I nd > < ns I ] (3.7) 


lns> and 1 nd> are IVannier state associated respectively 
with s and d band at the site n. 


g,s,d _ Center of the s,d band 

h^^ = Hopping term corresponding to s band 

(5 

hnm = Hopping term corresponding to d band 
T n ~ Hybridisation parameter for s-d interaction. 


In fact one can vi/rite down the matrix representation of 
in the follovj^ing forms 


H 






/ ^ 

Yi 

— 



1 ? 

~ I 



, «ds 

Hdd 

f 

1 

/ 


e? 

1 


/ 


H 


(3i8) 



The more explicit representation interms of t 2 g and 
symmetry components of the d— band is given by, 



H 1 
sd 


1 ‘^ds 

HpT 

«TE 

i 

i 

1 

! 

( 

HgT 

^EE 


It has been discussed by Ehrenreich and Schwartz (1976) that 
eventually all these would lead to a self-consistent CPA 
equation having the following matrix fo-rm: 

§ (z) = 1 + c(l-c) 6 fiz) [I-C( 1 - 2 c )6 + 

( £ - g(z))} F(z)]-^ . .... ( 3 . 10 ) 

where ^ ^ ~ ^ ^ B 

I = Unit matrix of (2x2) size. Here the energy parameters 
of the host B have matrix representation in s and d-bands . 
Explicit matrix representation for is shown in equation 
(3.8 ). 

The structure of the CPA equation (3.10) for the 
randan binary alloy corresponding to the Hamiltonian in equation 
( 3 .7 ) is same as the CPA equation in ( 3 , 6 ). Only difference 
is that here the self-energy S(z) and the disordered alloy 
Green Function ^(z) appear in the form of matrices having (2x2) 
representations in s and d-bands. Here in matrix form we have 
self-energy and Greai Function respectively given by the 
following expressions: 
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2 (z) = 


I S®(z) 


I 2 ( 7 > 


I 


i 


( 3 . 11 a) 




I 

i 




5 i z) 

t^sdr ' 


Sdd(-) ! 


(3.11b) 


v/here S®? S 2^^ are respectively the self-energies 
corresponding to the site-energies for s, d-bands and hybridi- 
sation parameter. It is very interesting to note that each 
of the Gre<3a Functions G^g(z)j andG^^(z) having 

representations in s,d-band indices is functional of all 
the component self-energies. Now in the k-space we have » 
hybridised s and d-bands model- Hamiltonian for the host, 
the crystalline system given by 


H( k) =. < k I H I k>= 


H ( k) 
ss^ 


di 


Ysd 


Hdd( k) 


(3,11c) 
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and G(k, z) = Spectral Function for the Crystalline System 


°ss 


’^ds 




G^ddc.z) 


where more 


explicitly in band structure terms 


have, for crystalline system 


G^^Ck.E) 


Vdd 


f ( ^ - V 33 S(S))( z -4-Vdd S(f)- 1ffdl 






I (z-e^ - ^'■^B'*^dd^^s^ 

'|)/''3S>-Ysd 1 


■ss ^yk), E) 


(3.12) 


where E (k) = e® -'-V^ S(k),V^^= hopping integral corres- 

S D ss 'SS 

ponding to s band for the crystalline system, 

= corresponding to d band for the same 
sd ~ hybridisation parameter 
S( k) = structure factor corresponding to the 
crystalline system. 
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Vbxxous Gxoon Functions fop t-ho cpystsllinG systorn « 


G (E) 


S (k, E) 

k 

+ 00 

’/ f is (1- E^(k)) dri 

-oo ^ 


" 1-00 
= " f 

-oo 

4'00 
= / 
-00 

2 E-i-lO 


z 


+ O0 

= 7 


-00 

E+iff^ 


%s (n, e)p^(ti) dti 

dd(^- 

- 4 - '"ss 'I- ^dd 

-t 

(’I- -Vsd] 

[ ( 2 -n) I Z - 0^ - (n- - Y sd i 

(3.13) 


v^e'oce' Pc('n) = density of states of s-band for the 

crystalline system, i.e, the host (B), 
Proceeding in the same way, one can have, integral representatior 
also for dr-band given by 


Sdd(^> 


-i-00 

-oo 


!z-eg -V 


(n- 4 ) 


ss 


V 


[{Z 


B 


V 




dd„ 


^ dn 


ss 


V 


dd 


}(Z -n) ] 


(3.14) 


where ^(3(”'^) = density of states of d-bend in the host 

similarly integral representation for sd Green Function 
is given by 
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G^dlE) = 


+ 00 

7 

-.oo 




[tz - e 


B 


V33 -SB)/VddHz-n) 


Y H ] 

’ sd 
(3.15) 


The relevant representations for disordered alloy Green 
Functions for various bands in CPA form are the following. 
Alloy spectral functions are given by, 


/ E-^o- VddS(k) , 


■1 


G(k, E) 


- S 


sd 


E- 2® “ V ^ S(k) 
o ss ^ ^ 


(3.16) 


More explicitly, 
Gdd^^»^^ = 


Ed(k)) 


[«E - S ® +(V,,/V^rt)i^n - HAk)?tE - +(E^-E/k) 


ss^'^dd^^ B 


0 ' B ‘‘r 


where E^( k) = B and structure for d-band 


sd 


] 


As before, 
Md* 


"b Md • 


(3.17) 


G^^(E) = / 




t'E -S^+(V33/VddX"B -ri )HE-S^+(ed_^-, 


(3J 8) 


Similarly integral representation for s band is given by 
; 1°° !E-i;^-tVdd/V3^r,-e|)}Ps(n).dTi 

g^^(E) = / — — 


-00 


[!E - -(V^yV3p(tl-e|)ItE SP - (ri-e 
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and 


G 



+C0 Sjd • 

[!E - 5;^ -( r,_ e|)!.E’-^ 

(3.20) 


We have generated the various Green's Functions in 

the integral representations by using smoothed version of d- 

band density of states for Cu through using terminator scheme 

by C.M.M. Nex in the continued fraction which we have discussed 

in detail in Chapter 2, Note that since s-band is featureless 

the standard practice to choose this is of semielliptic form 

as discussed by Hhrenreich and Schwartz (1976), We have 

chosen our disordered systems to be the substitutionally 

random binary alloy Cu^ selected the various 

parameters of Cu and Ni from the band structure Table (3,2a) 

et al, 

presented by A, Bansil/( 1 975) for Cu, Ni and the alloy 
composition, Cu^ see the d-band density of 

states for Cu, we would notice one thing that almost all 
the structures in the d-band is reproduced except the 
unusual height of the right most peak. This is due to 
finite size of the cluster and that we have to compromise so 
long as we are interested in studying the effect of local 
environment on the electronic states. The CPA equation for 
such a system in which s-d hybridisation is a important 
factor is basically a coupled non-linear equation and this 
can be solved iteratively starting from guessed solution for 



49 


K* ^sd as ( e| , -#), {Bq , ^ ) where 6' 

is a small complex part. The governing CPA equation (3.10) 

together with equations (3.18), (3.19) and (3.20) constitute 

d s 

the v.'orking equations to be solved for S ^ o' ^od' 
three self-energies and hence the three Green Functions 
G.dE), G (E), G -(E). So far as the numerical efficiency of 

qO $ S so 

S 

solving such equations is concerned, having attained 2^, 

Sq for a given energy E starting with guessed solutions, we 
use these as the initial iterative values at E + dE. At all 
energy intervals where these self-energies do not change 
violently, convergence is rapid. At several points convergence 
may be slow and upto 100 iterations may be necessary, 

3.2,2 Results and Discussions on CPA Density of States; | 

As stated earlier we have generated the d Green Function 
for the disordered alloy using the crystalline d- 

band density of states calculated by Recursion Method for a 
finite size cluster and ten level of recursion co-efficients I 
(a^, ) as shovn in Table (3.1), Haydock (1980) has ' 

discussed the Recursion algorithm in detail for the d-band 
of a f.c.c, crystal. We have made use of the most recent 
version of the Recursion Program Library in which the newly | 

I 

developed terminator scheme for Green Function by Haydock and | 
Hex (1984) has been used. The only mismatch is the right most I 
peak in the density of states in the pure ease compared to j 

I 

the density of states as shown by Bansil et al (1975). But | 



TABLE 3.1 


Chain P AkxMETERS for d Elecihons in an fee Crystal 



Orbital 1 Orbital 4 

0 00 0.100000£'+ 01 0.0 0.100000£+0! 

1 -0.736282£ - 02 0.295643£ - 02 - 0.262667£ - 02 0.211155£-02 

2 -0.285673£ - 01 0.235335£ - 02 - 0.833 1 32£ - 02 0.158428£-02 

3 -0.160718£ - 01 0.179840£- 02 - 0.962267£ - 02 0.121493£-02 

4 -O.I342I8£ - 01 0.162690£ - 02 - 0.540195£ - 02 0.601881£ - 03 

5 -0.112190£ - 01 0.156254£- 02 - 0.472927£ - 02 0.131840£-02 

6 -0.131124£- 01 O.174530£- 02 -0.398838£ - 02 0.869993£ - 03 

7 -0.67195g£ - 02 0.H3999£ - 02 -0.611713£ - 02 0.895845£ - 03 

8 -0.165180£ - 01 0.195674£- 02 - 0.21931 1£ - 02 0.854633£ - 03 

9 -0.123618£ - 01 0.173710£ - 02 -0.107605£ - 01 0.126003£ - 02 

•0 0 0 O.I68461£- 02 0.0 0.125162£ - 02 


Orbital 2 

0 • 0.0 0.100000£ + 01 

1 -0.7362g2£ - 02 0.295643£ - 02 

2 - 0.285674£ - 01 0.235335£ - 02 

3 - 0.203862£ - 01 0.205787£ - 02 

4 - 0.134712£ - 01 0.152931£ - 02 

5 0.907024£ - 03 0.16436]£ - 02 

6 - 0.596882£ - 02 0.155990£-02 

7 -0.463769£ - 02 0.13l672£-02 

8 -0.128520£ - 01 0.134450£ - 02 

9 0.]28076£- 02 0.1609I7£-02 

10 0.0 • 0.158227£ - 02 


Orbital 5 

0.0 0.100000£ + 01 

-0.262665£ - 02 0.21 1 155£ - 02 

- 0.833 1 3 1 £ - 02 0.1 58428£ - 02 

-0.1 16066£ - 01 0.121205£ - 02 

-0.12861 1£ - 02 0.610331£ - 03 

-0.736858£ - 02 0.135928£ - 02 

-0.194090£ - 02 0.985064£ - 03 

-0.532768£ - 02 0.829151£ - 03 

0.453 lg5£ - 04 0.955998£ - 03 

-0.123153£ - 01 0.127546£ - 02 

0.0 0.147460£ - 02 


Orbital 3 

0 0.0 0.100000£+01 

1 -0.736282£ - 02 0.295643£ - 02 

2 -0.285673£ - 01 0.235335£ - 02 

3 -0.203862£ - 01 0.205787£ - 02 

4 -0.134712£ - 01 0.152932£-02 

5 0.907012£- 03 0.164361£-02 

6 -0.596881£ - 02 0.155990£ - 02 

7 - 0.463773£ - 02 0. 1 3 1 672£ - 02 

8 - 0.128520£ - 01 0.134450£ - 02 

9 0.128075£ - 02 0.160917£ - 02 

10 0.0 0.158227£ - 02 


rest of the structure is similar, V'/e have repeated the CPA 
calculation for the hybridized s-d hamiltonian for various 
concentrations (c = ^ 23%, 40?^, 50%) of Ni in Cu hostC^ig"^'^ 

We observe as it is expected that the impurity sub— band is 
featureless in all cases and having almost Lorentzian shapes. 
As the concentration of Ni increases the tvvo major peaks in 
density of states of pure Cu start vanishing along with 
the emergence of featureless impurity band and finally for 
509^ Ni the impurity band density of states results in two 
big st3?uctur el ess humps. Note that the vanishing of the 
peaks in the density of states and the loss of structure in_^ 
impurity band is due to two reasons (i) artifact of single- 
site effective medium which tries to smooth out the local 
fluctuation effects (ii) effect of hybridisation. These 
results are in qualitative agreement v,/ith the old calculations 

. a in d 

by Stocks et al (1971), The self-energies j 2 ^0 

constitute the basic input for the calculation of the cluster 
CPA density of states. Our main goal vt’ould be to study the 
cluster effects coming from a pair of a tons coupled statisti- 
cally. The limitation of single-site CPA is recognised from 
real system results which constitute the basis for development 
of the generalisation of CPA and the new results that emerge 
out of it. 
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3.3 


Electronic Density of States of C-n Mi in Cluster 

. . 1— c "c 

CPA ' — — - — 

We have already given the introduction of ASF 

(Augmented Space Formlism) in Chapter 2 and also discussed 

the appropriate approximation vjhich leads to the self- 

consistent cluster CPA, At the initial stage ASF had 

proved to be successful in describing the scattering from 

clusters for single band model systems in 3-dimension by 

Kumar et al (1983), V,K, Srivastava (1982), Now, here, we 

are in a stage to apply it to a real system like Cu. Mi 

I c 

alloy whore one can judge the new structure in the impurity 
band that emerge from statistically coupled sites by making 
comparison with the CPA results that we have already presented 
in the previous section of this chapter. We have generated 
the various diagonal and off-diagonal Green Functions 
required for solving self-consistent CCPA equations in terms 
of Recursion Co-efficients suitable for f.c.c, d-band. 

The calculation we have pursued is a 2-site cluster 
CPA based on n general s-d model hamiltonian suitable for 
describing electronic states in transition noble metal alloys. 
The general tight binding form of the hamiltonian where s-d 
hybridisation effect is considered written in equation(3,7) 
in the previous section. We know that the d-band density 
of states for pure system has much more fine structures in 
noble and transition metals compared to s-band and it is also 
very sensitive to the hybridisation effects as discussed 



by Hodges et al (1966), So we give primary importance to 
the d-band of the constituents Cr and Ni of the disordered 
binary alloy Cu-]_^Ni^ with the realistic choice of centers 
of the tivo d resonances, i,e, the diagonal disorder of the 
d-band hamiltonian, off-diagonal disorder due to the vari- 
ation of the ividth of the d-bands and by generating appro- 
priate local Green Functions in Recursion Method for the 
d-band of a f,c,c, crystal. It has been seen as discussed 
in the previous section that d subband due to Ni impurity 
in Cu host is featureless in CPA and here we are going to 
investigate the appearance of structure over the impurity 
band for various concentrations (10?o, 23?o, 30 %, 40%, 50%) 
of Ni in Cu^_^Ni^ for 2-site cluster CPA calculations. 

Here in our calculation of 2-CPA density of states 
we have treated s-band and hybridisation in the ordinary 
CPA, While d— band we have treated in 2— CPA, We have alreacJy 
derived the 2-CPy\ results for a single band model Hamiltonian 
through equations (2,31 a, 2,31b, 2,35 and 2,36) in Chapter 2, 
Here we write down these main equations using the proper 
notations for various diagonal and off-diagonal Green Functions 
with respect to site induces corresponding to d— band. The 
augmented space graph which describes the scattering from 
statistically coupled with pair is an octagon as shown in 
Figure (2,2a^b) In our real calculations of density of 
states of Cu^^^Ni^ one should in fact come across the 
following self energies S®, E 2 ^ 2^®, 2^^, 2^^ ,xHere 



S z) 8 , 2 }^ 2 :)sre diagonal self-energies corresponding to the 
randomness in the site energy terms for s and d-bands 
respectively in the hamiltonian in equation(3 ,7 ) , Along 
with these we come across the two new off-diagonal self- 
energies 1 ^ gCz) 8.2^ ^(z) in cluster CPA that properly takes 
into account the randomness in hopping integrals corresponding 
to s and d-bands of the same s-d hamiltonian , As stated 
earlier we are interested in studying the diagonal and off- 
diagonal self-energies ^g( ,z) -(^hQ 2-CPA, The self- 

energies equations for d-band in 2-CPA has the following 
forms 




S S 




[ <^o^ 


o. ] h 


dd 

k* 


(3,21a) 


HereS^Cz) is the self-energy corresponding to the 
diagonal disorder in the Hamiltonian in equation (3 ,7) in 
its d-band term, 

V/e have here, j,k the site indices and 

H = Host d-band hamiltonian in augmented space 

<^dd 

h = Bond hamiltonian which connects the bond 
(of, If) with the rest of the octagon in 
Figure (2,2a, b) corre'spon-ding to d-t>and, 

g*^^z) ~ Green Function corresponding to the renorma- 

lised medium (6x6) obtained by partitioning 
the (2x2) matrix in the basis Iof>andl 1 f> 
from the (8x8) octagon, in d-band representation. 



(3.21b) 


S izY- 

1d^ 


H 


dd 

o1 


+ 


2 

a?^o 


2 


dd 


■*dd 

^oj '='jk 


^1k C ^ 


^dd 

*’ki 


where S ^ ^ is the self-energy corresponding to the off- 
diagonal disorder in the hamiltonian for its hopping term 
corresponding to d-band. 

Here, 

- Diagonal renormalisation parameter 

corresponding to the octagonal decoration 
Oj ( z) = Off-diagonal renormalisation parameter 
corresponding to the same octagonal 
decoration as discussed below„ 


and j,k being site indices as before. 

The renormalisation of the site and bond (in real 
space) of the rest of the lattice is carried out by the 
two renormalisation parameters o ^ and given by the 
following equations which we have derived in Chapter 2. 


"AA 


= °E£ = 


(2) 


meN^ “^km 


(3. 


AE 




Gw" (3 
km 


wheife, 


H 


( 2 ) 




and is the Green Function matrix corresponding to 

( 2\ 

H*" operator in the site indices k,m. 


22a) 
.22b) 



= Nearest neighbour of A 

Ng = Nearest neighbour of E 

P = Projection operator 

T = Transfer operator. 

The above equations (3,21a, 3,21b, 3.22a. and 3.22b) 

constitute self consistent equations for the self-energies S°(z) 

o 

and corresponding to d-band to be solved for 

characterising the translationally invariant effective medium 
in 2-CPA and hence to get the density of states. The 
appropriate Green Functions for d-band and s-bands with s-d 
hybridisation in which we have considered for calculational 
purposes d-band in 2-CPA, s-band and s-d hybridisation in 
CPA have the following analytic representations. 


( k,E) 


E-sjz)- S,U) 

#! 

-S3<j(z) S(k 


(3.23) 




,2CP 


'dd 


(k,£) = [. 


c 


o E^k)][(E- 2^ -2 ^ S(k))]-S 


2 

sd 


or finally, 


*2CP 

^dd 


(E) 


+ 00 

7 

-00 


Ce - 




- 


[!E-< +(V^3/V,X - >1 )KE 

2.^(P - £ P ) ■' 2 , 

B ^}-2^ 1 (3,24) 


'1d^ 


V 


sd' 


dd 



In the same manner, we have.s-band Green Function in 2CPA, 


G 


2CP 

ss 


(E) 


4 - 00 

/ 

--00 


{E-sgUSic/vJn- e|)iP^(n)dn 

[! E eUs = 

(11- e g)J - ^sd ^ 


(3.25) 





3.4 


Results and Discussion 


We have solved the CPA and the CCPA equations for 
Cu rich CuNi alloys in the concentration range ^ 0 %, 23^, 30% 

40% and 50% of Mi in Cu. Fig ,( 3.1- 3, 4) displays these results. 
Alloying broadens out most of the sharp features in the Cu 
d-band density of states, as do s-d hybridisation effects within 
CPA. With increasing concentration the impurity band grows 
at the expense of the states at the upper edge of the host band. 
Let us discuss the impurity band 10% Ni in Cu, against 
the Table 3,2b, which gives the eigenvalues of finite clusters 
of Cu andNi, At these low c aicentrations, larger Ni clusters 
are highly improbable. Within the CPA, details of features 
of clusters containing more than one Ni atom cannot be accurat- 
ely resolved. The featureless impurity band of the CPA may 
be associated with the eigenvalue around E=0,167 (row 1 of 
Table 3.2b) broadened by the immersion of the cluster in the 
infinite solid. The cluster-CPA begins to resolve features 
in the impurity band. The two peaks seen in the CCPA may be 
associated with the eigenvalues around E=0, 145-147 and 
0,189—190 (rows 2,3 of Table 3,2b) related to the bonding- 
antibonding levels (X3-X5) in Ni, The teirainology arises 
from the fact that in the tight-binding model the corres- 
ponding wave functions combine in the same way as bonding- 
antibonding wavefunctions in molecules. The small peak 
around E=0,06-0,07 probably arises because of the' energy 
level around those energies in small Ni clusters (rows, 2,3 
in Table 3.16 ). 
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State 

Cu 


Cuc.T6Ni«.,s(Vti“) , 

Cu<i.tNi,.s' 


— 

- r, 

0. 0161 

(0.0194, -0.0) 

(0.0194. 

-0. 0) 

(0.0226, -0.0) 

(0.0260, -0.0) 

0. 

cu{j:»' 

n 

0.S808 

(0.3822, -0.0038) 

(0.4408, 

- 0. 0050) 

(0.SS33, -0.0077) 

(0.3841, -0.0116) 


0.4510 

(0.4287, -0.0087) 

(0.4844, 

-G. 0140) 

(0.4112, -0.0124) 

(0.3966, -0.0144) 


j Fj!,/ 
hi ^ 

11 


(0.5645, -0.0286) 

<0. 5639, 

- 0. 0286) 

(0.5652, -0.0191) 

(0.5661, -0.0095) 

0. 


(0.5920, -0.0239) 

(0. 5950, 

-0.0197) 

(0.6150, -0.0147) 

(0.6349, -0.0070) 

Oj 

(Xi 

0.2178 

(0.2364, -0.0003) 

(0.2809, 

-0. 0002) 

(0.2587, -0.0007) 

(0.2884, -0.0014) 


Cu 1 

C. 2580 

(0.2765, -0.0006) 

(0.3290, 

- 0. 0005) 

(C.29S4, -0.0016) 

(0.3272. -0.0036) 


)x. 

0. 4997 

(0.4518, -0.0132) 

(0. 5003, 

- 0. 0206) 

(0.4236, -0.0150) 

(0,4022, -0.0155) 



0.5165 

(0.4581, -0.0148) 

(0.5039, 

- 0. 0222) 

(0.427U, -0.0158) 

(0.4037, -0.01 58) 


iX, 


(0.5317, -0.0310) 

(0. 5389, 

-0. 0327) 

(0.4969, -G. 0247) 

(0.4560, -0. 0188) 

0.3 



(0.5403, -0. 0313i 

(G. 5443, 

-0. 0327) 

(0.51S1, -0.0245) 

(0.4793, -0.0171) 

0.4 



(0. 6199, -0.0195.) 

(0.6303, 

-0. 0133) 

(0,6560, -0.0123) 

(0.6853, -0.0060; 

0,7 

Xi 


(0,6312, -0.0181) 

(0.6444, 

-0. 0117/ 

<0.6710, -0.0116) 

(0.7029, -O.OC58; 

0.7. 

• ^4- 

0.8372 

(0.B369, -0.0) 

(0.6369. 

-0. 0) 

(0.8366, -0.0) 

(0.836S, -0.0) 

C.& 

(i; 

Cu<Ll 

0.2347 

(0.2516. -0.0003) 

(0.2826, 

- 0. 0002} 

(0.2721, -0.0009) 

(0.2993, -0.0017) 


0. 3739 

(0.2780, -0.0035) 

(0.4361, 

- 0. 0044.J 

(0.3805, -0.0073) 

(0.3827, -0. 0113; 


(4 

G. 502 S 

(0.4531, -C.0135,J 

(0.5010, 

- 0. 0209) 

(C.424S, -0.0152) 

<0.4025, -0. 0156) 




(0.5293. -0.C3G4.1 

(0.5374. 

- 0. 03245 

<0.4939, -0.0240) 

(0.453&, -O.OiSS) 

0.3s 

Ni \Ll 


(0.5626, -e.0255i 

(0. 5621, 

-0. 0291; 

(0.5612, -0.0194) 

(0,56.00, -0. OD&S) 

0^ 5s 

4 


(C.6216, -0. 02 92, 

(0.6327, 

-0. OlSOi 

(C. 65B£. -0.0121) 

(C.65SS. -0.0060; 

0.71 


0. 6183 

<0.6187, -0.0) 

(0.6187, 

- 0. 0) 

(0.6161, -0.0) 

(0.6174, -0.0) 

0.61 


1.0016 

(1.0381, -0.0025) 

a. 0554, 

-0. 0016.^ 

(1.0722, -0.0031) 

(1.1045, -0.0021; 

1.33 

Ef'^ 

0.4086 

0.400S 

0.4582 


0. 394c 

C. 383:' 


£?■ 


C. 5700 

C. 5763 


0. Sf'O: 

i. 593 •: 

£.60 


(after Bansil et al 


1975 ) 


TABLE 3.2b 


( 

: Impurity levels found for a finite size cluster 
foimed by putting Cu atoms surrounding Ni 

atoms (starting from oneNi atom at the centre 
etc. ; 


3-2b 


No of Cu No of Nj Eigenvalue^Ryd) 


12 

1 

-0.0083,0.07, 0.167 

11 

2 

-0.0080,0.059,0.145,0.181 

10 

3 

-0.0080,0.049,0.145,0. 147,0.190 

9 

4 

-0. 0080, 0. 041, 0. 144, 0. 147, 0. 189 

8 

5 

-0.0083,0.033,0.144,0. 147,0.208 


1 


1 


-0.0031,0.12 
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Fig. 3.4 fd) Density of states for c=40V. Nl in Cui^fNic alloy in 2 CPA 
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As we go on increasing the Ni concentration, the 
structure develops strongly around 3C9o Ni in the impurity 
band. The structure is reminescent of the t2g structure in 
pure Ni and continues to be present upto 50^^ iMi, We may 
compare these results with the 13— size cluster calculation 
by Stocks using the Lloyd formula as reported by Faulkner 
(1982). They compare very well, except for the exaggeration 
of the rightmost " antibonding’' peak, which in our case is 
an artifact of the small cluster recursion used ( this is 
clear froui the similar exaggeration for pure Cu and Ni) , 

Such structure arising from clusters are not expected to be 
seen in the CPA. 

An alternative approach to the nearly free electron- 
tight binding (NFE-TB) type formulation we have followed in 
this work is the KKR-CPA developed for various metallic 
alloys (Stocks et al 1978j Temmerman et al 1978). The method 
has been made tractable by the development of fast Brillouin 
zone integration techniques (Bansil and Prasad 1975), In a 
review paper on alloys Faulkner (1982) compares the results 
on various binary alloys (CuNi, AgPd etc.) from the two 
approaches with experiment. His conclusion is that both 
methods yield as good a result on density of states. 

Gonis et al (1984) have attempted an embedded cluster 
approach (ECM) for AgPd, in which they embed clusters of 
different configurations in a KKR-CPA (single site) medium. 
The medium is not self-consistent with the cluster; Faulkner’ 
comments on the approach is instructive; "the method is 



just a cobbled up theory. It is hot self-consistent in the 
sense of the CPA or even the MCPA*’ , 

The marriage of thd recursion method with the augmented 
space techniques cut down considerably on computation times. 
Once the CPA has been generated, at most energies the CCPA 
requires about five iterations for convergence usinij rftulti- 
vatiate Nev;tcn— Raphson method. This is comparable to most 
KKR-CPA works reported. However at b#nd edges and where 
new structure appears in the CCPA often the convergence is Very 
slow requiring as many as fifty to a hundred iterations* Again, 
this is a common feature of all CPA works. 

In connection vath CCPA work often the statement has 
been made that' the fihe structure revealed in the cluster 
work is not ’’seen’* in say, photoemission experiments, whereas 
the simple CPA ( occassionally convoluted with a Lorentzi^ 
resolution function) seems to give a better fit to experiment. 
This fine structure was first discovered by Dean (196C,1'972) in 
computer simulation experiments. It was so startingly' 
different from the picture known at the time from crude mean 
field work that initially it was thought to be errors in 
computation’. However, careful analysis identified these 
structures as arising from lesonant scattering from clusters 
of different sizes and shapes immersed in the infinite matrix. 
The weight of the peaks are related to the probability of 
the clusters appearing in a random array, and is therefore 
strongly influenced by local short ranged order. The CPA being 



a single site mean field theory cannot reproduce such fine 
structure. Similarly experiments ( such as photoemission) 
always have finite resolving power. No optical experiments 
can measure density of states directly. Usually response 
involves the density of states convoluted v/ith functions 
(e,g, transfer matrices, surface and many-body effects) which 
for disordered systems may have strong energy dependences 
(Nowak and Dederichs 1985; Mookerjee 1986b), It is not 
surprising that the CPA, which smoothens out the density of 
states seems to agree better wtih experiments, which also 
measures the density of states smoothened out by different 
effects (Faulkner 1982). But this does not mean that cluster 
resonance structure is unimportant. Indeed, the electronic 
driving force that would cause an order-dis order transition 
in a crystal has its origin in the resonance structure in 
the density of states. 



3*5 Conclusion 

^ - ■ 


It would be very interesting to note that CCPA in 
augmented space formalism has been also applied to describe 
phonon density of states in disordered binary alloys like 
NiPt (1985), NiPd and. NiCr alloys (1982) where the traditional 
mass CPA failed because of strong clustering and force 
constant, disorders. The underlying physics of cluster 
induced states in density of states of electron and phonon 
have close similarities in the sense that wherever the 
diagonal and off-diagonal disorder or the off-diagonal one 
is strong one can do calculation in CCPA to describe the 
density of states within ASF which is an unified approach 
for random systems. 

In conclusion, we have presented a method of obtain- 
ing a self-consistent cluster CPA. Here, we have for the 
first time performed the calculation of electrical density 
of states in CCPA for realistic disordered metallic alloy 
like Cu^_^Ni^, a canonical system which has received attention 
both from theoretical and experimental point of view in the 
past. Our calculation is done in 2-site cluster CPA and 
resembles the new structures in the density of states which 
appears to be of very general nature. 

This shows how important is the cluster effects at 
its minimum size and in fact for a system where the disorder 
is strong one should do CCPA where the size of the cluster 



concerned should be of the order of the number of nearest 
neighbours, Ofcourse, the formalism is valid there but 
the numerical evaluation of the various self-energies would 
be very cumbersome. 



CmPTER 4 


RESIDU/\L CDHD UCTIVITY OF m P OM B IMA RY ALLOYS 
iiT trac CC P}\ “■ 

4 • 1 _Intrqcb-.-P-'to^,Y .. Qsjqa jrks on tho DovoIop raGnt of Rasi dual 
C oHd uc^b ivit y ^ 

Wg have concontrated so far on tho oloctronic deiosity 
of states in disordered alloys within the co-herent Potential 
Approximation (CFA) and its generalisation CCPA in which 
effects of statistically coupled clusters as well as off- 
diagonal disorder have been treated self-consis tently. 

Another important physical quantity of interest from 
theoretical as well as experimental point of view is the 
residual conductivity describing the electronic states and 
its transport in disordered metallic alloys in d,c. electric 
field. Our interest would be to study the conductivity as 
an enerc'y functional throughout the band in the spirit of 
CCPA in ASF, The residual conduct-ivity would be the value 
of this functional at E = E^, the Fermi energy. 

We shall first recapitulate the older theories of 
conductivity based on the CPA in detail with a view to 
point out the shortcomings like the vanishing of the 
vertex corri,3Ction, non- random nature of the current etc, 
which can lead to significant loss of fine structure in 
the conductivity results as compared to the CCPA results 
based on ASF, Such a recapitulation will indicate a v/ay 
of generalising the older theories so as to over come the 
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shortcomings and formulate a method that is feasible and 
yields phys ically valid results in realistic systems, 

4 . 2 Various Dev. 5 lopm en tal Stages for the Res i dua 1 C on du c t iv i ty 
in Disordered Arioys ; ” ^ 

Theoretical development of conductivity in random 
alloys has been carried out in several simple situations. 

The Ziman-Faber theory (Ziman 1961, Bradley 1962, Faber 
and Ziman 1965) developed originally for liquids but 
applicable also to alloys, is valid for weak scattering of 
the conduction electrons by a distribution of dilute random 
scatters rs. The Ziman-Faber theory a nearly free electron 
picture within the Bom approximation is valid only for low 
concentrations of scatterers, the resulting mean free path 
is much larger than Debroglie wavelength there. For the case 
of strong scattering, in transition or noble metal alloys, 
one may use the extended Ziman theory which replaces the 
scattering potentials by equivalent t matrices is discussed 
by Evans et al (1971)., Roth and Singh (1982) pointed out 
that the approach neglects multiple scattering effects. 

The principal multiple scattering formalisms are the 
average T-matrix approximation (ATA) and the CPA, the former 
beirig nonself-consis ten t whiJ-e the latter incorporates self- 
consistency, It has been discussed by various authors 
Ducastelle (1972, 1974), Muller-Hartmann, E, (1973) and 
Mookerjee (1973) that CPA maintains the Herglotz analytic 
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properties of the Green Function for all concentrations, 
energies and scattering strengths. Since CPA deals with 
single site scattering it does not take into account the 
back scattering effects which is not small in strongly 
disordered systems. 

The calculation of electrical conductivity for dis- 
ordered binary alloys in the spirit of CPA was introduced 
by Velicky et al (1969) starting from Kubo-Greenwood 
formalism (1 957,58) .The CPA formulation of conductivity has 
been also applied to s-d model Hamiltonian for simple model 
systems having semi elliptic density of states both for s and 
d-bands. But in both the two cases only diagonal disorder 
was taken into account. Later Niizeki (1974, 1977) tried 
to generalise the formulation for both diagonal and off- 
diagonal disorder but in the spiri't of CPA and was able to 
show that the vertex correction is not zero even in a 
single band case. Until recently Butler et al (1984) have 
made calculation of electrical resistivity of Ag. -Pd alloys 
in the spirit of KKR-CPA theory, a CPA theory derived from 
the multiple scattering methods having many features in 
common with the KKR (Korringa, Kohn and Rostoker) approach 
to band theory. The basic KKR-CPA equations were derived 
by Seven (1970), and later on Shiba (1971) and Gyorffy 
(1972) put them in convenient forms which are of general 
use today. The general formulation on electrical 
conductivity based on the Kubo-Greenwood formalism is also 
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carried out by Butler (1985) through application in one 

i 

dimension. 

Our (1984 1986) aim is to present a general formu- 

lation of residual conductivity based on Kubo-Greenwood 
form of the conductivity expression for three dimensional 
m.odel systems and a. realistic system like Cu^_^Ni^ in the 
spirit of cluster co-herent Potential Approximation (CCPA) 
via the ASF, It has been established that ASF is not only 
proved to yield analytic results for density of states of 
disordered alloys but this is also true for the residual 
conductivity in CCPA of the same,. 


4,3 Older CPA Theories of Residual Conductivity 

4 3-, 1 Theory of Static Electrical Conductivity for Randcro 

Bjna ry iO)j?o vs "by VeTrc^ 

The; theory of static electrical conductivity for 
substitutionally disordered random alloys has been carried 
out by Velicky (1969) starting from Kubo-Greenwood form 
of the conductivity expression. It is based on the single- 
site, CI^A suitable for a tight binding form of the Hamiltonian 
as in equation (3,1a ) in Chapter 3, The starting point 
of the conductivity formulation is the Kubo-Greenwood 
expression ( 1 957 , 1 958) 


<cr 


.oc^ 


m Q ' 


pP 6 (t]- H) > 


(4.1) 
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where 


< > 


P 


a,p 


= Notation used for configuration averaging 
= Component of the linear momentum in a given 
direction ( or P ) 


m 


D 

f(ri,T) 

H 


Effective mass of the electron 

Volume of the unit cell 

Fermi-Dirac distribution at temperature T 

Tight binding Hamiltonian of the disordered alloys. 


Mow , 


< k 1 p'^ i k' > = 


m ( — -) e ( k) 6 1 

9 k® kk 


m va ( k) 

kk’ 


(4^2.) 


Now using the operator identity, 

- (27tir''[G( TT) - ] 

± + + 
p * -► _Ti — iO 


(4.3) 


Velicky expressed conductivity expression in a convenient 
form for calculation' in the CPA, in the following manner, 

= 2|ot iP“lf 

with = (4^)'"^ Tr p®[K(T]''", p^.,!!"") + 

“ K(h"*, pP, h'*’) - P^fTi'^) - K(t]“, I? ,p“) ] 

i 4- 

T) •* ± (4^4b) 



where 


(4.4c) 


K( , p, Z 2 ) = < G( z^ )p G( Z 2 ) > 

is the vertex function, G( z^ ) being the Resolvent corres- 
ponding to the random alloy Hamiltonian, sVe note that the 
quantity K ( z.] , p, Z 2 ) to be averaged is a product of two 
random operators Vi/hich can not be factored as 
<G(z^)> p <G(z 2 )^ in general. 

Simplification of the equation (4,4c) in terms of 
configuration averaged Green Function can be done by the 
following manner. The standard equation defining T-matrix 


is 



G 

= G + G T G 

(4.5) 

where, G = 

configuration averaged resolvent 


T = 

G = 

functional of self-energy 2(E) 
'unaveraged Green Function. 



Averaging equation (4,5) on both sides, <T> = 0 for the 

self-consistent 2(E), Now substituting equation (4,5) in 
the expression for vertex function K in equation (4,4b) 
v/e have, 

K = < G p G > 

= GpS + S<T5p5T>G 

G (p +r) G (4,6a) 

where vertex correction is defined as, 


r = 


< T 5 p G T > 


(4,6b) 



Velicky was able to express the equation for K involving 
vertex correction T in a convenient form suitable for 
calculation in CPA and he showed that vertex correction 
vanishes for a single band Hamiltonian having short ranged 
potentials in the following way. 

The expression for K is v^ritten as, 


K = Gp'G + G 2^fn><nIKfn><nlG 

(4.7a) 

using r = S ^ I n > < nl K I n > < n I (4,7b) 

where cv is a function of ( whose explicit form is 
shown by Velicky as, 




< \( ) t^( > 

[ 1 + F( Z-, ) < t^( z^ ) t^( Z2)>F( Z 2 ) ^ 


(4,8a) 


where, scattering matrix t ( z) 


= [s n ^ -2 ( z)] F( z) 5 


(4,8b) 


Alloy Green Function F(z) 

— < n I G( z) 1 n > 

= F'’ (z -S (z)) (4.8c) 

Here, 

= site energy; ^ ( z) = site diagonal self energy 


F°(z) = Green Function corresponding to the pure system. 



Now taking diagonal element of equation (4,7a) with 
respect to Wannier basis, one can write, 


< m I K I m > = < ,m ) G p G I m > + 

LVn 'Kin > 


(4.9a) 


where «( 2 ., ) 

m— n 1 


< m I G( z) I n > 


<m I gH z-2( z))m>; GH z) = \ lie >fe-e ( k)]“^ < k I 

ix 


(4,9b) 


Velicky obtained the solution of equation (4,7a) by 
expressing R.H.S, in a Fourier Transformed form as, 


< m I K 1 m > - N' 


where m ® 


■1 \ 

k ® 
-ik,.I^ 


a. 

k k 


-1 


< m 1 G p G I m > 


-ik.R 


and Ajt = 2 e 

K m 


■m 




(4.10a) 

(4.10b) 
(4.1 Oc) 


Now, <m 1 G p G f m > 


1 

N 


It: 


-• 13 ( Z.J ) - e ( k)] ^ m ( k) 

(k)]-^ 


[ Z2 - S( Z 2 ^ “ ^ 


(4,11) 


By time reversal symmetry. 



Therefore, we have, from equation (4,11) 


<mlGp5lm>=0 (4.1 

Hence, from equations (4,10a, 4,10b) and (4,7b) 
we have, = 0, 

The vanishing of the vertex correction is the 
artifact of single— site CPA theory for short ranged 
potentials, 

Velicky derived the CPA equations for residual 
conductivity in a very simple form after proving the 
vanishing of vertex correction in the following manner. 

The quantity vt/ithin the energy integral in equation 

(4,4a) can be evaluated by simplification of the expression 
for vertex function K in equation (4,6a), We start with 
the simplification of the term 

K(ri^, , r)”) 

= G (n'") p^ a(Ti-) 

= s Ik _ e(k) - 2 (ti)] '''' < k I 
k 

pP S i k’ > [^'^ -e(k) -S(T 1 '■)]"■■* <k’l 

k' 

== S S !k>(G^k, E) + i gV, E)) 

k k’ k k’ 

(G^ (k* , E) - i (k‘, E))<kq 

= Slk>{p^ [G^^ (ic, E) + G?^ 'k,E)]3<k I (4*1 
k 

where G^ = Real part of spectral function G( k, E), 



= Imaginary part of spectral function G(k, E) and 
Til = £ + 

In the same manner-, * 

2 2 

f?, ^ ^ 1 k>^pP G.^ ( ic, E) ■+ G^ ( k, E) 3 } < k I- 
k 


(4.13b) 


^ ' 2 2 

=5:lk> { p^Cg'"^ (k,E) - G^ (k,E)l><kI 

k 


and 


(4.13c) 


K {11 , p'^.tl-) 


-Z -♦ 


= 5 Ik ■->CpP[G^ (k,E) - G^" (k,E).] } ■ <k I 

K 


(4.13d) 


Therefore., finally, 
a j3 


• p p 


( ^ ., h) 

a 3 ^ 2 ^ 

■2^ Tr [ 2 1 ic > 4 G-^ ( k, E) < k J ] 

47ir 

rvin ■* **l 

[ ^n ^k ^ 


7t' 


„a„P j2 - 

C ^i; (k, E)l. 


7C 


k 


(4.14a) 


Note that <n f k> and <k In > is simplified by using 


r k>- N“''/^ Sm 


I m > 


(4,14b) 


where |k> = Block-ketj I m> = Vifannier-^ket, 


Now, finally the conductivity expression has the 
following form. 


= 2®! T 6ri{-M) 5 v“(t) /(k) [lmG(k,^)]2 

,10-00 Sn ^ 



Velicky was able to derive the weak scattering and 
Boltzman limit from the above equation (4.15) in straight- 
forward manner, 

■S£§J:^^ §..'L..Qonductivity for a s-d Model Hamiltonian, 

cM by Srouers and Ved^^y^ 

In the previous article we have discussed how the 
static electrical conductivity of a single-band model with 
short ranged random scatterers was shown to be exactly 
solvable in the CPA. Brouers and Vedyayeb (1972) made an 
attempt to extend the CPA theory to the calculation of conducti- 
vity of a two-band model introduced by Levin and Ehrenreich 
Levin et al,(l968). Levin ( 1970 ) 

(19701 f /( lAhich neglects the widths of the d— levels due 
to d-d hopping, produces alloy density of states exhibiting 
hybridisation gaps at the d~levels. The previous theories 
by Friedel (1962) and Gomes (1966) which takes into account 
the effect of s-d scattering on the resistivity of transition 
metal alloys was considered only for dilute alloys, 

Brouers and Vedyayeb (1972) were successful for the 
generalisation which takes into account the, width -the d-band 
( of— course within a semielliptic density of states model) 
and the effect of s-d hybridisation so that it can be used 
for the complete range of concentrations'. It provides a 
useful tool for investigating possible deviation from 
Nordheim c(1-c) law (c being concentration of the impurity) 
behaviour of the resistivity of transition and noble metal 



alloys as well as the influence of a-d hybridisation. On 
the other hand calculation of conductivity by Brouer .-S. Vedyayerv 
(1972) is able to reproduce the weak scattering limit 
and dilute concentration limit as derived by Gomes (1966) 
for transition metal alloys, 

Brouers et al tried to calculate the configuration 
alloy averaged electrical conductivity of a disordered binary 
AcBi_c whose Hamiltonian has the form, 


„alloy 

n 


I E^{k) I Imdxndl 

m^ 


+ 2 
n 


where, 


n 


^s,d 

BZ 


mn 


i k R 

f nd><nd i+ S y («^)( I k ><nd i e “ 
k BZ ^ ® 


(4,16a) 


= 

\ = 
(k) = 


Site energy corresponding to d-band at the 
nth lattice position 

Lattice coordinate corresponding to mth site 


o 


and-s true ture corresponding to s-band or d-band 


= Brillouin zone and ^ ^ 

= N-1 5 E'^fk) 


(4,16b) 


Now, as it is expected conductivity will be a matrix 
in the { 1 I } representation . -they o- obtained 

the s-s, s-d and dd representation of the conductivity — in 
the following manners 



^ dri(- Tr < V Im Q,(i1+ iO) ^ Im G( tt+ iO) 


(4,17) 


Total conductivity, 


^ + <^ , + <T , + cr 

ss dd sd da 


(4,18) 


Note that the operators % and g of the above equations 


are (2x2) matrices in the basis of s and d bands, . Brouers 
et al derived the follov/ing results for s-s, s-^d and d-d 
components of conductivity , ' j 

or,, = -2^1 / dn (-11) I v^(k) [ini'G -dV T 1 + iO)]' 
ss TCQ . * ^ dr]-' k S ss ’ 


(4,19a) 


0 j. dT(-||) I v^ (1^) [ Im G^^(k, r]+ 


(4.19b) 


^ds = 2 . f^/dnC-l^)' |d^Ck) ;^(k) 


x[lm (k,Ti+ iO)]^ (4.19c) 

In the present case, ■fchd'''S;^,lf-energy has to be 

defined in the representation of d states whereas Gj^ can be 

i ' ■ . . , 

written as , a (2x2) matrix in the basis t ^ ^ 


z - Eg(k) 


.■.;,G(k;z) = 


■f ,.-V -Y" 


Z -. e^(i) -“ 2^1 


(4,20a) 



The three Green Functions, z), G^^(k, z) and G^^(k, z) 

are as follovi/s: 


Gss(k»z) = C z - E^( k) -y^^Ez -2 ^ ( 1 ) ] 

( 4 . 20 b) 

G^^(k,z) = tz - Sq - A~i^) - y"^[ 2 -2 =(k)]-'' 5 

( 4 . 20 c) 

and 

Gs^(k, z) = Y" { [z - Eg( k) ][ z - 2 ^ _ e k) -Y 

( 4 . 20 d) 

where, v = Drift velocity of electron in the s-band 

aE^dJ) 
a ?■ 

e Yk) = ae ®( k) and = Drift velocity of electron in 

d-band 

3k' 


a = constant proportionality factor. 


As in the usual CPA equation one can have the self- 
consistent condition for d-band 


with 


<t^> - 


t^ = re^ - r 1 (e^.2^^F (-2^)] 

n n O'' n o' -^dd ^ o'^-' 


dM-1 


( 4 . 21 a) 

( 4 . 21 b) 


The final CPA equation of the self-enerqy 2 ^ is. 



- Sq) (z, Sq) (eg - 2^) (4,22a) 


2 ?^) = 


where e 



+ (1-c) 



(4,22b) 


and, (2,2^) = Q{2'^r^ S(Pk iz -S^-e"^ ( k) - 

(4,22c) 


Y = hybridisation factor , 

Equation (4,22a) for self-energy 2^( z) is solved 
iteratively with the choice that in the pure metal both 
s and d unhybridised bands have the same shape and are 
related by a scaling factor , The results of density of 
states and the position of the Fermi level for various 
concentrations have been shown in Figure (4,1^b) , The 
corresponding calculation of conductivity is also shown in 
Figure (4^1 c). 


Results and Discussion on the Residual Conductivity Calculation 
of Brouors e t , 

As stated earlier that the effect of hyoridisation and 
the finite width of the d-band is taken into account to overcome 
limitations of the orevious theories by Levin et al (1968,1970); 
Levin (1970); Levin an d Ehronreich (1971) on static residual 
conductivity. Although density of states that has been used for 
d and s~band of the pure system is assumed to be of semi- 
elliptic form yet the results show up the effect of hybridi- 
sation and finite width of the d-band on residtial resistivity 
for noble/transition metal alloys atleast qualitatively. 

If we focus our attention on resistivity 



(P) versus concentrations (c) curve in Figure (4, 2a^ ) on®-' can 
see the striking deviations from' the linear dependence of 
resistivity versus concentration for small concentrations 
and from Nordheim’s c(1-c) law for hi^er concentrations. 

In Figure (4,2a) (resistivity vs, concentration) the region 
of low concentration is magnified in the upper part of the 
figure (4 ,2b ) and shows the change of the slope of the 
resistivity curve, correlated with the minimum in the density 
of states. The asymmetry of the total curve and the sharp 
maximum at ' c = 75^ emphasizes the significant contribution 
of s-d hybridisation in that concentration range, 

Mott (N,F, Atott and H, Jones — 1958) has emphasized 
that in transition metals and transition metal alloys the 
resistance is mainly due to scattering process in which an 
electron makes a transition from the s to the d— band, the 
probability of such a transition being proportional to the 
density of states in the d-band. With this picture Mott 
has explained qualitatively the striking deviation from 
Nordheim’s Rule (Pa c(1-c)) exhibiting a Matterhorn type 
behaviour for rather large concentrations (50 to 70%), In 
Figure ( 4 »1c) conductivity versus E curve the three contri- 
butions to conductivity are plotted. Here, 

one can clearly see a very striking feature that c , and 
can be larger than in the majority band regime even for 
a concentration c = 15%, This shows the strong influence 
of s— d hybridisation in the conductivity. 
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Fig. 4.1 (a 8, b ) The d density of states of the two s-d alloy model 
for various increasing concentrations 
(c) Totol dc electrical conductivity at T=0 , o^c csd 
and conductivities versus Ep ’ 

(from Brouers and Vedyayev ,1972 ) 




Total residual resistivity versus concentration in units of 
its maxirmum value. The asynrtrr^etry of the total curve and 
the sharp maximum at 75 7# emphasizes the Importance 
of s-d hy-bridizatlon In that concentration ranae 


The above results that emerge from Brouers Vedyayev 
calculation are in agreement with Mott’s prediction and 
experimental results by Nellis et al (1970) which shows an 
abrupt change in slope of resistivity at about lb% of U in Pd, 

4.3,3 Theory of Electrical Conductivity for a Disordered 

Hamiltonian having “bobh Diagonal and Off-Diagonal 

Disorder ' — — 

We have so far discussed conductivity for an Hamiltonian 
involving diagonal disorder only in single - site CPA, The 
effect of off-diagonal disorder in a model tight binding 
Hamiltonian for the electrical conductivity was originally 
investigated by Niizeki and Endo (1974), They showed that 
the conductivity is very sensitive to the off-diagonal 
disorder as well as the diagonal disorder.. On the other 
hand, Fukuyama et al (1974) investigated the case where 
transfer integrals are additively related. Several authors 
also tried to generalize further to the case of the model 
tight-binding Hamiltonian introduced by Blackman et al 
(1971) when no restriction is imposed on the transfer 
integrals, A remarkable property of the Hamiltonian is that, 
in the dilute limit, only the s-like partial waves (those 
with the complete symmetry of the point cfroup of the lattice) 
of the scattering due to an isolated, impurity have non- zero 
scattering amplitudes, and, hence, the total scattering 
amplitude is isotropic. Consequently, the vertex correction 
as well as the backward scattering integral in the linearised 



where, I^p = Unit matrix operator. Green Function corres- 
ponding to the random Hamiltonian is given by 

G(n,E) = (G(z) - H(n) + H( n, z) (4,24b) 

Now, G( z) = <G(n,z) > 

= x .G^^ (n,z) + yG®(n,z) (4,24c) 

where x and y denote respectively concentration of A and B 
components in the A-B alloy. An alternative form for self- 
consistency, given by 

< t (n, z) > 

= x t'^(n,z) + y t®(n, 2 ) = 0 (4,25a) 


where^ t(n,z) is the single-site T— matrix given by, 

t( n, z) = { - [H(n) - H(h,z)]G(z)?“'^ [H(n) - H(n,z) 1 

(4.25b) 

The relation between the t(n,z) and the random Green Function 
G(n,z) is given by, 

G(n,z) = S(z) + G(z) t(n,.z) G(z) (4.25c) 


Niizeki defined the current operator I and the position 
operator r^p in the following way; 


r 


op 


S ! ■ n > Rp < n I 


(4,26) 



I = 


(4.27a) 


where ICn) = " 

^ where e = electric charge. 


■(4.27b) 


Now, the static electrical, ;.con due tivity, given by Kubo- 
Greenwo’dd fomuia,:' ,! ) , ; 


cy:^ -= 

ap 0 


Tr- < (^I — H) I 6 (z i — H ) > 

a ^ op : - B ■ op 


(4.28a) 


where, a ' ,"p"'are spatial direction indices,- Q = - the volume 
of the unit cell and z is the' energy with ah infinitesimal 
imaginary part O , Conductivity expression can be ' ' 

<• ■■'' (- 1 ) ( 2 )' 

decomposed -Into two parts it .and ~ b-; a§ :foIlew,si-, ^ 


ba'= - 'b -h b(2 

ap ap bP 


(4.28b) 


Here, 


Tr <Ig G^(z) Ip Gjz) > 


(.4.28c) 


Re<Tr Ip G^( z) > 


(4,28d) 


Where G^(z) corresponds to Grosai Functions for z+ iO « 
Making use of the identity. 


G^(z) I G^(2) = 

( 2) 

a „ can be written as 
ap 

„( 2 ) ^ e% 

‘^a p IxQ 


K 'h 


(4.29a) 


Re Tr <(m;p:j^p G^(z) > 


(4,29b) 



1 


(4.29c) 


where. 


^"■opIb 


3 ^ H( k ) 

3k Sko 
a P 


Effecti v e Cu rrent Operator : The Verte.v Function is defined as, 

K(a} = G,(z) I Gjz) (4.30) 

Now, Niizeki, introduced an idea of effective (or coherent) 
current oporcator J(^) which is periodic and energy dependent 
but yields the same vertex function as I, i.e,. 


K(z) = <G^(z ) J( z) g _(2 ) > (4.31) 

which means <G_j_(z) [l- J(z)Jg (Z)>= o (4.32) 

Taking transpose one can also write, 

<G_(z) [ I - J(z) ]g^.(z) > =0 (4.33) 

using. J(z) = and inserting the following expressions^ 

I = J(z) -p [I-J(z)] 

- J(z) + [I _ J(z) ], 


Novy, the conductivity comes out of tho following form, 
('i) ^ ^nc 


o 


ap 


ap 


ap 


( 4 ,34a ) 


Vitiore, 



= 4- Tr[J^2)<G_^(z)Jp (ii) G_{*) >] 

= ir Tr [Jj;a)' Kp( Bj ] 


(4.34b) 


and 


inc 

ap 




- G+(s)[ Ip - Jp(z)]G(z) > 

(4,34c) 


c ir^ o 

According to Niizeki ( 1 977,) cr«:p , respect-iyaly. -ooherent 

and incoherent parts of the conductivity. 


It is important to note that J(2) defined by equations 
(4,31) and ( i,32) can be evaluated easily in the CPA prior to 
evaluating the vertex function in the following way as 
discussed by Niizeki and Hoshino (197?) , 


Neglecting cluster effects on the current term as before 
in CPA for self-energy, Niizeki introduced the effective 
single-site current operator. The underlying philosophy 
remains same and self-consistency is imposed on single site 
current operatorj(n ,z ) which is the analog of the 

single-site energy. 

We have, 

G(^) = ^ J(n,z) (4,35) 

The effective single^site current operator J( n,e ) is 
determined by the condition 

<G_^(n,z) [ l(n) - J(n,z) ] G_ (n,z)> = 0 (4,36) 

using equation (4,32). Now using expression for G from 



equation (4,24b) and writing the zeroth site for J(n,z) 


J(0,z) = < [l^p - 1 - t^(0,z) G^(z)3 i(o) 

[1op G (^) tJO,z)]x“<t^ (0,2) 

G_^(z) J(0,z) G_(z) t_^ (0,z) > (4.37) 

where t_^ (0,z) is the Oth-site t matrix for energy have 
been infinitesimal imaginary part o"^, Now, the equation 
(4,37) provides the evaluation of J(2 ) in the CPA to the 
solution of finite system of linear equations for the matrix 
elements of J(0,z ) between the basis vectors. 

Evaluation of Vertex Correction 

The vertex correction T (a ) in the vertex function 
K(z) is defined by the equation 

K(z) = (z ) [J(z) + r (2 ) ] G - (z) (4,38a) 

Velicky (1969) has shovi/n that the vertex correction 
can be divided in the CPy\ into single-site vertex corrections 
as 

r(z) = 2 r(n,z) (4.38b) 

Niizeki and Hoshino (1977) showed that the single-site 
vertex corrections satisfy a coupled set of equations given 
by 

r(n,z) = < t^ (n,e) Gjz)(J(2) + r(n*,z)) 

G_(z) t_(n, 2) > (4.39a) 



Equation (4,39a) for the representative site 0 can be 
written as, 

r(0,z) = < t^. (0,z) K(Z) t^ (0,2) > 

- <t^.(0,z) G^.(2)r(0,2) G_^(z) tj0,2) > 

(4.39b) 

The above equation gives the representation of vertex- 
correction in single— site terms. The vertex function and the 
vertex corrections are determined through introducing an 
effective current operator J (z) different from J(^ ) by the 
equation 

K(2) = Gjz) J*(z) G_(z) (4.40a) 

Now, 

j"" (a) = S /(n,2) (4.40b) 

and 

J (n,z) = J(n,2) + r(n,2) (4,40c) 

Now, from equation (4,39a) adding J(n,z ) on both sides and 
after' some algebraic manipulation one can get the single- 
site representation of the effective current operator 
J (n,z) . 

J (n,2) as shown by Niizeki and Hoshino (1977 ) is 


given by 



J*(n,z) = <[ 'op + V 

lCn) .[l„p -i- G_(^) tjn.z)] > 

+ <t+Cn,z) G_^(z) (p?^ An'? ) ) a-(==) t_(n,z)> 
/ (4,41 ) 

and the vertex function K(z) given by, 

K^( 2 ) = - < G_j_( n,2 ),[ Ji ( 2) + I( n) ® J ( n-, z )] G__(n, z) .> 

■p':'/'-.: ■; :'vP, '' ^'(4.42) 

N ow , eQua ti on ( 4 ,41 ) and (4, 42), with ecjua ti •n (4,40a,b,c). can 
be regarded as a self-consistency condition for J (n,z). 

So, the above mathematical outline shows that in ... 
general vertex correction can not be zero in CPA and in fact 
it can be evaluated in the spirit of effebtive current approach 
as stated above. The brief account of numerical results for 
conductivity calculations done by Niizeki and Hoshino (1977) 
in a model system is given here, 

• ■ ' . ' V ' 

i... ^ J ’■ ' 

Numerical Investigation of Static Conductivity by Niizeki 
an^ ITo sHirTo — - 

We have previously discussed the theory of electrical 
conductivity of a random binary alloy which is described by 
an one electron Hamiltonian in the tight binding approximation' 
with additive randomness. Such type of Hamiltonian can be 
divided into the unperturbed Hamiltonian describing a periodic 
reference medium and the scattering potential which can be 



divided .further into single— site potentials. The main 
difference between the previous CPA theory of Seven (1967) 
and the CPA theory used by Niizeki and Hoshino (1977); • 

is that the single-site p'otentials are long range and nonlocal 
in general in the latter case. The long range potential may 
give rise to anisotropy in the differential cross section 
and since may lead to nonvanishing of the vertex correction. 

Previous numerical investigations of the electrical 
conductivity in the CPA of a random alloy have been carried 
out only for the case where the scattering potentials are 
short ranged (Velicky 1969, Levin et al 1970) or the case 
v;here the long range can be reduced to the former case by a 
suitable transformation (Niizeki and Endo 1974j Niizeki and 
Hoshino, 1 975 ) , The vertex correction vanishes in these cases, 

Niizeki and Hoshino (1 977)n . calculated numerically the 
effects of the vertex correction and the nonlocality of the 
potential on the conductivity. 

The model they considered is described by a tight 
binding Hamiltonian given as, 

H = S h (li > <j I + f j > <i T ) 

< ij > 

where the summation is limited to nearest neighbour pair 
of sites in a simple cubic lattice, I i > indicates an 
atomic orbital of the ith site and h^j is the transfer integral 
a random variable which takes care of the two possible 



X Y 

negative values h and h with possibilities X and Y 

respectively. This model has only off-diagonal disorder 

but the character of the off-diagonal disorder is different 

from that in the Blackman-Esterling-Berk model ( 1 971 ),B:la'clcman 

(1975). The present Hamiltonian corresponds to the bond 

model in the classical percolation problem and the latter 

one to the site model, Niizeki and Hoshino (1977) computed the 

density of states and conductivity as an energy functional 

X Y 

for X = Y = 0,5, h = -1 ,0 and h = -0,5, Computation of 

conductivity cy and its three components cr and 

as energy functional for the above alloy parameters chosen^- 

shew very significant effects of vertex correction,! as shown in 
Fig .4 .3a, b) . 

They also computed the total conductivity o 
without vertex correction. All these results are shown 
in Figure (4 ,3ajl^ , They pointed out from their numerical 
results that the conductivity is decreased by the vertex 
correction, because the differential cross-section of the 
scattering of an electron due to X or Y type bond embedded 
in the periodic medium in the CPA is larger in the backward 
scattering (the initial and the final wave vectors k and 
k’ are antipaigallel) than the forward sca^terjjjig (k and ^ 
k' are parallel) , This is eviden in the weak-scattering 
limit when the scattering amplitude is given in the Born 
approximation is proportional to, say, Cos (k + k’ )a/2 

whose absolute value is larger at k = -k’ than at 

-♦ -♦ 

k' = k. 






Drawbacks of the Old CPA Theories _on Static Conductivity 
and the Subsequent Improvements ;' 

In the previous section- (4*3 ) we have given the 
mathematical outline of the old theories on static conducti- 
vity, Drawbacks are the following; 

(i) Single-site CPA theoi'y of static conductivity 
by Velicky proves the vanishing of vertex correction (i.e,, 
the back-scattering effects) for short- ranged potentials of’ 
the alloy components. This assumption is not true in 
general as the potential may be long ranged, 

(ii) Brouers & Vedyayeb (1972) generalised the treatment 
Velicky for a s-d model Hamiltonian, But they performed 

the conductivity calculation for a very crude model system 

vdiere they assumed the s and d bands to have semielliptic 

electronic density of states. They pointed out only the 

effect of hybridisation in the residual resistivity results 

(i,e,, the deviation -from Nordheim' s Rule), 

(I977a,b,c) 

(iii) Niizeki /first proved that the vertex connection 
in general is not zero in CPA, His calculation was for a 
model simple cubic lattice. But his conductivity results 
give a good deal of information about the dependence of 
conductivity on vertex correction and randomness in current. 



4,3,4 S tatic Electrical Conductivity of Random Binary 

Al loys in A SF s 

Recapitulations of the old theories in electrical 
conductivity provided us the idea of studidng the effect 
of vertex correction and, the randomness in the current on 
the conductivity for general Hamiltonian where off-diagonal 
disorder as well as hybridisation should be taken into 
account. Our aim would be mainly to study the effect of 
vertex correction and the randomness in current on static 
conductivity in CCPA via ASF, The effect of off— diagonal 
disorder and statistical clustering are incorporated through 
the self-energies corresponding to a finite size of the 
cluster within CCPA through simplification of the configuration 
averaged conductivity expression ( Kabo-Greenwood form 1957,58), 

4.3.5. 

An Augmented Space FprmAlisit. , (.ASF ) of Electrical. .Cond u ctivity 
in Random Alloys wi t h in a Clust er Co -herent P otenti^ 
Api^ roximat ion (CCPA}T~" ' ' 

The aim of this section is to develop an explicit 
formalism for calculs. ting the conductivity in the CCPA within 
the augmented space formalism, v7e shall illustrate the 
methodology by application to a simple three-dimensional 
model alloy first and finally to a realistic system like 

Just as the Herglotz property is essential for any 
physically meaningful approximation for the one-particle 
propagator so that the reality of the eigen energies and 



positivity and single-valuedness of the density of the states 
is preserved, Mookerjee (1976) showed that for any approxi- 
mation for response functions, the vertex correction must be 
related to the self-energy via a Ward identity, so that 
averaged macroscopic conservation laws are not violated. 

Given the self-energy in a CCPA, the .7ard identity then 
provides the vertex correction. This is closely related to 
the variational procedure of Popeilawski (1979) and that 
followed by Roth and Singh (1982), The vertex correction 
with the self-energy is sufficient input for the Bethe— 
Salpeter equations for the response functions. The Dyson 
and Bethe-Salpeter equations for the CPA were obtained by 
laborious diagrammatic resummation by Leath (1970). However, 
within the augmented space picture, in which k-labelled 
states are scattered by configuration fluctuations, both 
were obtained in the CCPA byMookeriee (1975b,c), These 
form the background of our present development, 

We have already stated that there have been earlier 
attempts at the theory of conductivity within which off- 
diagonal disorder was taken into account (Niizeki 1977a, b,c, 
Koshino and V/ataba 1977, Miizeki and Hoshino 1977). There 
are several essential differences between the treatment of 
these and allied works and that of ours (Thakur et al, 1986), 
All these earlier works v\/ere within the framework of the 
single site CPA, gehieralising, for long-range potentials, 
the pioneering work of Velicky (1969). Single-site CPA is 



intrinsically unsatisfactory for problems with off-diagonal 
disorder, where at least the 2 CP/., should be considered as 
the lowest-order approximation. Moreover, except for cases 
vvhe.re is either the arithmetic or the geometric average 

of V,A and the CPA equations themselves do not simplify.. 

In fact, all CPa v;ork '-‘/ith off-diagonal disorder simply 
ass’jmes this to be true. As we shall see later in realistic 
alloys, the presence of short-range o::der invalidates 
these assumptions, VJithin the augmented space approach, no 
such assumption is necessary (Kumar et al 1982). There is 
nothing in these works which indicate how to generalise to 
the CCPA without encountering analytic dif fic ul 'ties , The 
treatment of the effects due to randomness in the current 
terms is also different in the two approaches. Finally, 
all these earlier works envisage only ofi'-diagonal disorder, 
appropriate to the problem of bond or site percolation? 
whereas, our approach is applied to a model of an alloy 
in which there is both diagonal and off-diagcnal disorder, 
statistically coupled. It would be instructive to see 
how to deal v'ith such a model within the earlier formalisms. 


The .Mathematics of Stati c C onduc ti vity .Formulation in ASF ; 

We sliall take as our starting pointan one-electron, 
tight binding Hamiltonian with one orbital per site 


where P^^ 


H 


2. e.P. 

1 11 


i > <i f and 



?V(r, -r.)T,. (4.4 3) 

I i><j I , The basis functions 



are a complete set of linearly in<5eoendent site labelled 
functions. In general, orthogonality is unnecessary, 
although for conputational pu^rposes it is convenient, A 
useful basis are the pseudo-ordoi tals of the chemical pseudo^ 
potential approach of Bullett (1975), as used by Mookerjee 
et al (1933;. Extension to situations with several orbitals 
per site is sti^aight forward, involving replacement of 
scalar quantities like e and V by matrices [e^^], etc, 

where n , m label the orbitals , 7e may .go beyond the 
tight binding and extend our calculations to muffin-tin 
type of calculations. These extensions will be discussed 
towards the end of this section. 

In a completely random system we have 


e(rj_) = 


e..,n. 
A i 


!(1 


Hi) 


jj 


'A 


'B 




V(r, 


T' ' 


n 


(4.44) 


= Vgg + V.| (n. -I- n.) -1- V n n 
wh-.ra. V, = + Vg3 - ZV.gj ^ . 


Mere n. is a random occupation variable for the site r. 

J- 

by the species A, ’.7ith probability. 


Pr(nj_) = c6(nj^ - 1) -i- (1 - c)6(nj_) 

where c is the concentration of A in B. 

■Viookerjee (1973) showed that the configuration average 
of any function of { n^} may be expressed as 



F({n^}) = <flF({MJ)If> 


(4.^; 


'.'/here F is the same functional of the operator as F 

is a function of the n^ , hh is related to the probability 

density of n^^ . In our case, it is an operator in a vector 

space. of rank two, ’-vhose representation in a basis 
, i i , . 

{Vq, V,) is ■, 


M = i I • 

“1 lc( 1 -c) lV 2 i_e J 

^ N 

F is an operator in the space # = 9.^ of rank 2 and 

ip . ^ ^ 

lf>= ^ 'reference’ configuration, a member 

in the basis of $ , 

This formalism has been described in great detail in 
several publications, already referred to, and has been 
utilised to generate Ilerglotz, self-consistent CCPA averaged 
to the single particle propagator G( r- , r.|E) via the self- 
energy matrix (Kumar at al 1982), 


The starting point of the conductivity calculations 
is the Kubo expression 

0-“^ = dE(-df/dE) <y^\E) ^.47) 

—do 


'.vhere f(E) is the Fermi function and 

(E) = 2'rt-h Tr{j“^(EI - M)/<3(EI _ H)} 

= (■h/7i:)Re Tr{/ G(E''’)/g(E'‘) -j°‘G( E^) / G( E"^) 1 
= Cs^^E'^E”) - e'^)}, (4.48) 



G( z) is the single particle propagator or the Green operator 

cc cc 

and j =2 is the current operator. Like the off- 

diagonal part of the Harailtonian , vve may write 




a 


a 




a 


Hi n . 




'■NX th 


; jyCR) = 3^(R)+ 3 bB^ R)-2j^B^ H) 


and 


j^(R) = (ie/h)R ( R) or ( k) = (e/h)<^^ • 


( 


In a random system vve have to configuration average 
o’( E) and this we shall do within the augmented space 
formalism. Both the operators G and j are random, involving 
the random variables Cn^^} , In the earlier attempts by 

Velicky ( 1969), and Brouers and Vedyayev (1972), since no off- 
diagonal disorder vms taken into account, randomness in 
j did not arise. In the vvo 3 rk of the Japanese group, random- 
ness in j was accounted for by finding an effective energy 
dependent current j(B), much in the spirit of the coherent 
potential itself. The approach of howak and Diederichs (1982), 
in the can text of the analogous case of neutron scattering, 
is also of this nature, Fo3.1owing equation ^45) we obtain 

<S°^^( , z^) = (fi/^)^e Tr<flj G(z^)J^ G( z^) I f>; G=( zI-H)"""^ 


with 


(4.50) 



H = Hg I-K3a-^b> 1 ‘>1 « '''1 IJ ''-'iTiJ ®'' 

(4.51) 




and 
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j ^bbAj « I ^ 

ro 

M. = 

I 

I \X . , , 5b- , , 

(N> 

II 

I 

\5c;> I . 4 « % * , 


The expressions in (452) raay be regrouped slightly* 
The advantage of such a regrouping would become clear 
immediately. If vjq write Ih == cP^^i m^^, then 


where 


a 


a 


a 


a 

.a 

Jb 


if 3 'hhj X'Ff + 3b T.J (mi mj )+ j^ij ®m.j ) 

(4.53) 

o cc , p a , . ct 

('’“^) J BB 2c(1-c)j^g 


a , V a , V ct 

oJaa + (l-2cJjAB - 


BB 


. cc 
^c 


2. . o: ^ .a 
c (a 


M ■' JbB ^ 


On ) 

•Q / ■ 


(4.54) 


ct a. a 

j is the averaged current, while j, and j are corrections 

due to the disorder in the current terms, arising, out of 

off-diagonal disorder. The last of these terms is closely 

related to the short-range order in the system. If we form 

an alloy by quenching the melt from a high temperature T^,. 

then the configuration which is at equilibrium at T^ is 

frozen-in at the lov/est temperature. Follow3.ng Clapp and 

Moss (1966, 1968), the short-range order parameter, 



a ( R) = / 


d^k 


D e' 


-ik,r 


8 




1 -i-2c(1-c)V^(k)/kn,T 


D = fd^RC^(R) (4.55) 


and V^(R) = -i- V^g - unless V^( R) = 0, 

o: ( R) 0, fcr xR ^ 0 and there is always short-range order 
in the alloy. The Clapp and i.ioss treatment clearly indicates 
that in the presence of off-diagonal disorder, unless the 
alloy is quenched from a temperature much above th.e critical 
order/ disorder transition temperature where the short-range 
order vanishes, v»’e must have shor't- range order characterised 
by ( R) . Our subsequent treatment which includes off-diagonal 
disorder but ignores short-range order is therefore, not 
st.rictly consistent. In any realistic calculation for real 
alloys we must incorporate short-range order? Kaplan and 
Gray (1977) have indicated how to do this within the augmented 
space f ornalism. 


Incorqo orating the partitioning as in equation (4,5 4) 
into (45)0) we obtain nine terras from the cross terms between 
and Consider a few of these; 


ap 

S-j (2^,22) 


(.t 


a 


,-P /•- / 


xRe Tr<fIj3G(z^)f^G(z2)!f > 




. .'x3( r^ - r^) If > 


Re / (dh/8n^)j^Ck) l ' 


P(1) / 




( 4.5 6 ) 



aP , . _ 


-P, 


(h/Tt) Re Tr<f G( 7.^2) 2 z^ ) I f > 

= (h/’t) Re [c(1-c)]V2 j“{V?k) 

X <fiIG(Vv"^1^Ja^Gj,-r„)G{r„-?^;J2)If > 


= (h/K) Rs /(d2k/8n3)[c(1_c)] j“(k)r^2)(z, .ZySk) 

( 4 . 57 ) 

'’/here \ f^y = I v^> '.x vj > . 

Sg^Cz-jjZ^) = (h/Tt) Re Tr<f Ij^ G( ) j^G( 2^)1 f > 


= (h/it) Re G“(ri-rp 


X <f IG(Jrt;zi)jy;rt)G(r„-r. jzj,)l fi> 


= (h/7t) Re J'(d®k/8i>^)[c(1-c)]'/^ 

( 4 . 58 ) 

The first term is the familiar one which arises because 

a, . 

of the averaged current It is completely diagonal 

a 

in the augmented space. The next two terms are off- 
diagonal in augmented space, being taken between the two , 
different configuration bases I f > and • These nine 

contributions may be grouped into categories; those 
diagonal in augmented space and those which are off-diagonal. 
Let us now examine these contributions in some detail. 

The Bethe-Salpeter equation for the two particle 


propagator 



( 4 . 59 ) 


G (r^. 


r. 


“k^ 1 


) = [gC 




:,JG( r. , r. 




)]. 


may be schematically expressed as in Figure (4,5al The 
equation for L' ^ ^ (z-|»Z 2 l^) maybe obtained by contracting 
with ( r^-r^ ) shom schematically in Figure (4.5b) . 
The resulting Bethe-Salpeter equation may be written as 
shovm; 


*^ ' ''( z., , Z 2 ? k) = G( 5 k) j^( k)G( Z 2 ; k) G( z-j j k)G( Z 2 ;k) 


X / ( d^k'/8'n?)A( z^ , Z2J k, k' )L^ ^ ^ \ z^ , Z2; k' ) , 

(4.6 0) 


This equation provides a self-consistent calculation for 
Lb '^(z^,Z 2 ?k) provided the vertex correction a( z^ , z^ 5 k, k' ) 
can be found. 

The solution to the vertex correction follows from 
the Vi/ard identity 


6 S( z^ , Z2? k) = T ( d^kV8Tc^)/-( ,, Z2 J k, k' ) 5G( , Z 2 ? k’ ) , 

( 4.6 1 ) 


Let us examine the structure of the vertex correction 
term in the CCFA scheme in augm anted space,. Let us focus 
our attention on the sca^'-tering diagrams in augmented space 
for the vertex correction. These diagrams were discussed in 
detail in J.lookerjee (1975b,c) and we shall take up our 
analysis from that work. 



The single-site CPA involves only those diagrams 
which include one site (Figure 4.4a ), so that the vertex 
correction is totally diagonal in real space and, Jcecause 
of homogeneity in the augmented space, is k indepaident* 
This leads to the vanishing contribution to conductivity, 
arising out of the antis^imimetry of j (k) for reflections, 

3 . , 

around the origin (Velicky 1969). In the two or any larger 
cluster CCPA, the vertex term is neither diagonal in real 
space (Figure 4,4b ) nor in k— space. However, inspection 
of the diagrams reveals an important topological property. 
Any diagram may be decomposed into tivo parts (Figure 4,4c) 
such that either part is any of the diagrams of the series 
shown in Figure ( 4 .4d.) , Algebraically this implies 


A ( 


zigjk, k' ) 


X ( z>| ;k) K 


Z 2 I k' 


( 4 .62) 


The separability of the vertex kernel allows us to solve 
the integral equation (4,61) explicitly. 





(E jk') 


(4.63) 


A(E'^,E'*'5k,k' ) 
^■] (Ef k) 


= ^2*^ ^ » k) ^2^ ) (4.64) 

= ImE(E.ik) 

{/ (d^k'/STt^) [ira G(k’ ,E)][lm2(E,k')YT^ 

a'(E;g) 

{j (d®k'/8’^)E'(E;k' )[2'(E;lc' )-1 .olG^CEjk') 


wl'i ere 21 ' 


Energy derivative of 21(E5k) 


( 466 ) 



Let us* concentrate on a tivo-site CCPAj the Bethe- 
Salpeter equation (4,60) may be rewritten as shovm in 
Figure (4,5 q ) ,On the right-hand side are throe types of 
contributions s 

(1) The first term is of zeroth order and its 
contribution in k- space nay be written as 


G( z. 


-»• 6 -> 

k)j3(k)( 


7 ® Ic 'z 


A( z. 



(4.67) 


(2) The second term ai-ises out of sing'le-site 

scattering. Note that in the single-site CPAj, this is 

/ 

the only vertex correction tenr that remains. In the 
2CPA tooj this term must be taken into account. Note, that 
its contribution in the 2CPA is not identical to that in 
the ICPa, since all the internal propagators are G(z,k) 
calculated within the 2CPA and not the 1CPA. However, the 
structure of the terms is the same in any CP approximation. 
The vertex is totally diagonal in real space and independent 
of k. Its contribution is 

/L^^^^\z^,Z2|k)(d^kV87r3) (4,68) 


( 3) Finally there are the terms arising out of 
two-site scattering. These are the lowest-order terms 
of the socalled 'maximally crossed' diagrams. These 
terms are totally absent in the single-site CP.b with 
short-range potentials. The contribution of these terras 
in real space is given bys 



S E 
s, s ' 


G( z, ; r^-r^ )G( Zg; r^, -r^ M ( z, , Zjj 

\ z, , Zj;r^,-r^) (4.69) 


Taking the Fourier transform, we obtain the contribution in 
k-space, 


A( z. 




, Zr,| k) // 


d ^kV d^k» 
8^ 87C^ 


;k-i-k'-k»)X(z2”,k")L^^^\z-j,Z2' ■ 


(4.70) 


If v/e look at equation (4.63-1*66) , we note thatS(zjk) = 2^ E^s(k) 
where s( k) = S exp(ik,x) summed over the nearest-neighbour 
vectors x on the lattice. The vertex terms arising out of 
these diagrams contain a part which is k independent (arising 
out of E^) and a part that does depend on k. We shall 
absorb the k-independen t term with the second term (as it 
has the same structure). The remaining vertex correction 
term thus arises essentially from the off-diagonal part of 
the self-energy. Thus, 


, 6 r Ad^k’ d^ic: 


Si ( Zi , Z 2 ) = f j^( k) jg( k)A( , ZpJ 

0 TC^ 


(er^)^ 


jg( k)A( z^ , Z 2 J k) A( z., ? k-i-kf-k'’)A( Z 25 k" )L^ ^ \ z.| , ^ 


(4.71 ) 


Equations (1 -(1 ,71) provide us with the basis for obtaining 

the conductivity, arising out of the averaged current. 


t 



-■ / S 


q 7 i 3 


JAO;'-!' Z' ' 375 ^ [A'i’ , M-' V (3 

^l^vhz^.z^Sk) = ■"nal6'Wr^^1-^25‘'> 

r^X.lrW') 


d^k' 


^T]a, a6^^1 * ^2^ 6a • 

(4.72) 


For computational purposes, v/e may project the matrices on 
to vectors and supermatrices into matrices. Take, for 
example, a case with two orbitals per site. Then the 2x2 
matrices of the kind may be written as the 4-vector 


6gp = (6G^^,6G^2»^21'^^22^ 


T 


( 4 . 73 ) 


vvhile suoermatrices like A ^ may be written as 

f a6 


A„ 


i ^11 , 1 1 

I ' 

CM 

*» 

A 

11 ,21 

A ; 

11,22 1 


2,12 

■'^ 2,21 

A i 

12,22 I 

1 “21 , 1 1 

^21 , 1 2 

^21 , 21 

■^ 21,22 

I "'^ 22,11 

J 

l||» 

^' 22,12 

^ 22,21 

^' 22,22 1 

and ( 4 ^ 1 ) 

may now be 

rewritten 

as 


L ^ \ , Z2Jk) -:- 2 ^pq^ 2 :^ » ^2 ' 

cjr 

X n d^k787i?)Aq^( , z^^k, k» )lP^ ^\z^, z^J k* ) 


(4 .7 5) 


I 



and 


6sp(z^,Z2;k) = 2 / (c[^kV8TV^)Ap2(z^,Z2Jk, k')5g^(z^,Z2|k>) 

(4.7)6) 

Equations (4,75) and (4,76) novi/ foim the basis of our conductivity 
calcula ti ons , 


T ki. ,Wga.^.^iS i t ttri,n.g Limit 

Before we contemplate application of our formalism 
to any system, let us examine the weak scattering limit and 
compare our results mth the well Icnown results which 
already exist in this regime. 

Let us first express the current j_(k) as j(k) 

a 

cos'j'( k, i) where j(k) is e3E(k)/9k and9(k,,r) is the angle 
between the direction of k and a unit vector in the direction 
of r. Moreover, we shall write 




» ^2 


k) = 


- .p 


k)A( z 


1 » 


k)'.V( z^ ,Z2;k) (4 77) 


so that the Bethe-Salpeter equation becomes 


W(z^,Z2|k) = 1 -I- /( d^kV87f^)A( 
X cos 9 ( k, k’ )W( k’ , 


z-j » Zp I k, k' )a( 

2i , Z 2 ) 


Z2?k‘ )( j( k’ )/j( k) ) 
(4.78) 


For z^ = E*”, Zp = E” we may identify 6 E(e'^, E“; k) = -2Im2(E,k) 
-1 1 

= rQ(4sJ"'an isotropic life time. From the Ward Identity we get 
'rQ^(k) = -2 '^(d^kV8^)A(E'^E“jk,k' )lniG(E,kM 


(4.79) 



‘ Noting that 


G(E''',k)G(E";k) =_2ir~'(k) ImG(E;k) 


(4,80) 


and in the weak scattering limit the self-energies are small, 
so that Im G( S, k) is strongly peaked around the energy 
shell E = E( k) . The dominant contribution to the integral 
on the left-hand side of equation (4,78) then comes from this 
energy shell. Using these approjdmations we get V/(E',E”jk) = 
'1' ( k)/'i^Q( k) where the identify a diffusion life time 



= -2/(d k'/87i^)A(E^,E''5k,k» )lnG(E,’k' ) [1-cos 


(k,k')h 
(4.81 ) 


This is the Rubio result that reduces to the Ziman formula 
(Rubio 1969). 


Turning to the Kubo formula, if we neglect the 
corrections due to randomness in’ the current (which is small 
in the regime under consideration) 


XX/ N 

<r (n) 


(e%)/-^ j''(k)[j^k)lG(E,k)l^ /V(E‘^,E";k) 




x( 1 ). 




-Re ‘ ''(E- ,E''4k)3. 


(4.82) 


Because ff directional isotropy = c^^and 

ja(k) - 3( k) , Again the dominant contribution to the 

integral comes from the energy shell H = E( k) '. In this 
shell the main contribution comes from the first term within 



\ I / 


the integral (Rubio 1969). _Velicky (1969) points out that 

i 

the neglect of the last teim is essentially valid ini the 
regime where the Boltzmann equation is valid itself. This 
approximation should thus yield the Boltzmann result. 
Evaluating the first teira on the energy shell and noting 
that we are in a temperature regime v/here kgT is Small 

3 

o = (2h^-/3)/ [-Ini3(E,k)]T(k)j^(k)(_df/dE)] 

QTtr 

= ne^(Ep)/m , (4.83) 

This is the well known Drude result. 

Comparing equations (4,79) and (4,81) we notice that the 
vertex corrections lead to the inclusion of back-scattering 
processes.- Velicky (1969) and Brouers and Vedyayev (1972) 
have shovm that since the effective self-energy in the single- 
site CPP. is short-ranged and the off-diagonal part only 
nearest-neighbour, the back-scattering vanishes and there is 
no contribution from the vertex corrections. When off- 
diagonal disorder is present and we go beyond the single site 
CPA, this is no longer true and we have to discuss the 
vertex corrections, 

Applicgi.ti..on .to .^a .Simple Model Case : 

Before we go on to apply our formalism to realistic 
models of alloys, let us illustrate the procedure through 
an application to a simple, yet non- trivial model case. We 



I J o 


shall consider only one orbital per site. We must consider 
a three-dimensional lattice, as in any lower dimensions all 
electronic states are supposed to be localised and no 
cluster CPA can be vsxpocted to reflect this. Vie shall choose 
a simple cubic lattice. Our choice of this particular 
lattice is entirely for ease of calculation. The conductivity 
computation involves several Brillouin zone integrations* and 
the Brillouin zone for a simple cubic lattice is also simple. 
This makes the integrations stjoaightforivard; The input 
Green functions are calculated by the recursion method. 


We shall choosea 50-50 alloy, so that the corrections 
because of the randomness in the current (being functions of 
c(1-c)) are maximum. Other alloy pararaeters are 
®A - Vgg = 0.33; V^/Vg 3 = 0.5; V^g/Vgg = 0.75, 


Figure ( 4 , 6 <?) shows the conductivity v/ithin the CCPA 
as a function of E . The, dotted curve is the conductivity 
with vertex and random current corrections, while the full 
curve shows the uncorrected conductivity. It is instructive 
to compare this v-/ith the numerical results of Hiizelci and 
Hoshino (1977;. This also involved a 50-50 alloy on the 
simple cubic lattice, and with • 0.5, but it 

involved only off-diagonal disorder, with e„ = e™ = 0. 
Firstly, because of our different e^ and e^., the conductivity 
is not s\mmetric around E = 0 like that of Niizeki and 
Hoshino ( 1977a, b,c). The uncorrected conductivity roughly 
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, -‘t, • K ^oiise of the randomness 

The corjrections which arise becau^ 

, , . . . , entirely analogous 

in the current involve equations wnicn ai'- 

to (4.6Q-('1-.71) , but v/hich involve self-ene^^D^®® S(z,k) ihich 

are of f-dic.gonal in the augraented space, language 

^ J .1 to • j as corrections to 

of graph tnoory, these self-energies ari^® 


nted space 'decorations*. 


the host Hamiltonian Hg because of augraenx^^ 

Let us take as an example, the 2CPA. The diagonal element of 
the self-energy in augraented space was partitioning 

the 8x8 cluster matrix (corresponding to S - 2x2“ configurations 
of tv/o sites and their four possible q on orations a.-^, AB, 


BA and DB) 


into that spanned by Of and If thau spanned by 


be obtained by 


the rest. The off-diagonal elements may u, 

partitioning into a part spanned by Of and Of^ (for example) 
ivhich will yield < f fE^^( z) I f q> j or by Of ^nd 1 f which will 
yield <f E^(z)IfQ> and so on. Having these self- 

energies of f-diagOnal in augmented space? '''^® proceed with 
them and equations (4^0) to (4,71) exactly n® oefore. 

It only remains for us to indicate bow the procedure 
is modified if we have more than one orbif®^ site. It 
is impo.i'tant to do this, as in any realisf^® situation we 
will be involved, with more than one orbital site. The 
Ward identity and the Beths-Saloetor equafl®^® now look as 


follows: 
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reflects the features of the density of states;. Figure 
shows the correction due to randomness in the current. 

Like the corresponding term ©f Wiizeki and Hoshino 

( 1977 ), it is of varying sign over the range of Ep, 

Figure (4 sho'vs the vertex correction contributions.. 

It is always negative, the reason for which, based on the 
Boltzmann transport theory, has been discussed in some 
length by the above authors (Niizeki and Hoshino - 1977) 

These corrections are appreciable in some energy ranges. 

Its maximum is the centre of the averaged band En = 6 = e 
Again, the reason for expecting this has been discussed in 
detail by Niizeki and Hoshino (1977 ), At its largest, the . 
vertex correction is about 20^4, an estimate in agreement 
with the work of Niizeki and Hoshino (1977 ), 

It should be noted hare that although the vertex 
corrections lower the conductivity, it would hot be 
feasible to look for signs of localisation in small cluster 
CPA* The reason for this is clear if we return to the 
analysis of atleast one whole- class of diagrams, namely 
the maximally crossed diagrams in augmented space. The 
various CCPA only include the first few terms of this infinite 
class, and cannot therefore reasonably reproduce any 
information about localisation. However, in those alloy 
systems v/hera v/e do not expect localised states near the 
Fermi level, our formalism is a useful one. 




(d) 


Fig. 4.4 (a) Scattering diagrams in au^ented space for 

the vertex coming' from single - site scatterir^. 

Cb) Two -site scattering. 

(c) Separability ef the vertex diagram s- 

(d) Diagrams for tie function X { 2 ; k }. 
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g. 4-5 Ca) EHagrammatic representation of the Bethe - 

Salpeter equotion for the two particle propagator G 

fill} n 

(t>) Derivati<Mi of L ^ from G by contraction 

The diagrammatic r^resentation of t»e Bethe - 
Salpeter equation for 






{c} Correctiorjs irt ccmductivity arising out of 
the randomness in current 



We apply our formalism in the next section to more 

realistic models of sd trans iti on/noble metal alloys, 

incorporating vertex corrections and using realistic^ 

densities of states generated by the recursion method of 

Haydock et al (1972), In several such alloys, like 

Cu Ni. 1 theband widths of the constituents are quite 
C I c 

different, necessitating the introduction of of f-c-iagonal 
disorder. This is intended to supplement the initial 
qualitative fomulation of Brouers and Vedyayev (1972), 

4.3,6 y\n^ Augmented Space Formulation of ,th e Electrical 

Co'nductivrtv for CiL iT^i ' thin ’a“tirus'’ter Coherent 
' . t-c c 

Pote r^ar Alopro^dmation 

In the previous section (Mookerjee at al 1985) 

the theory of static residual conductivity of substitutional 

disordered systems is discussed in detail and application 

is made to a single band model on a simple cubic lattice 

in the spirit of cluster coherent Potential Approximation 

(CCPA) (Kumar et al 1982). It has been shown there. that 

the randomness in the local current and the vertex correction 

may play important roles in conductivity. It has been 

discussed by various authors (Kiizeki ~ 1977a, b,c) 

Chitnavis and Leath (1983) 

Chitna’vis and Paghavan - 1983^that the vertex correction 
has a non-vanishing contribution to the conductivity not 
only for systems with off-diagonal disorder but as well 
as for multiband systems even in CPA. So the vertex 



correction for realistic ipodel systems should be treated 

with care in discussing electronic transport properties. 

In this section our aim is to discuss the static residual 

conductivity for a s-d transiti on/nob le metal alloy like 

Cua Ni in the light of the cluster coherent Potential 
1 —c c 

Approxima ti on ( CCPA ) . 

Vi/a have already discussed that calculation of 
static electrical conductivity for s-d system is cai-ried 
out by Brouers and Vedyayev (1972) to a very simple 
model system where s and d bands are represented by serai- 
elliptic density of states. Their calculation v/as carried 
out in CPA and importance of vertex correction of off- 
diagonal disorder is ignored. Most of the earlier attempts 
made in static conductivity is not well suited for real 
systems because of the adlioc model for the system chosen 
or for the representation of the density of states by a 
crude analytic form. On the other hand we have chos®i 
Cu^_^Nic as a realistic system- and the various band para- 
meters-, are extracted from reliable band structure 
calculations (Bansil et al - 1975). Our aim is partly 
to illustrate the usefulness of the augmented space method 
for realistic calculation for given parameters of Cu and 
Ni; V/e do not expect spectacular effects in the CCPA 
for this. Mere, the vertex correction shows up to such 
an amount that changes of the conductivity results from 
that of without vertex correction being not much, yet it 



can not be ignored. Havever, it is a, convenient starting 
point for development of the augmented space techniques 
to conductivity. The importance of vertex corarection for 
low and high concentrations is discussed in detail through 
a generalised treatment suitable for multiband system:. 

In a previous section (Thaku:: et al 1986) vve 
performed the calculation of electronic density of states 
in the augmented space formalism for Cu^_^Ni^ for various 
concentrations , I7e have treated the scattering from a 
pair of sites having both diagonal and off-diagonal 
disorder j, The two self-energies (diagonal) and 
( off-diagcnal) obtained therein are the basic inputs of 
our present calculation, V\fe are able to get fine structure 
and effect of vertex correction both at lov/ and high concen- 
tration regimes through the green functions having a 
functional dependence on self-energies. 

The most interesting feature that emerges here in 

conductivity results as an energy functional is that 

conductivity resembles most-ly the 2CPA structure of the 

density of states. The difference of c endue tivity results 

betvveen witli and without vertex correction in Cu^ Ni 

I -c c 

is the appearance of a dip in the middle of the impurity 
band regime in the foxmer while the latter re serabl'es mostly 
the structure of 2CPA density of states. 



Multiband Foimulation for Conductivity in Augmented 
Spaces 

The suitable Hamiltonian for Cu-Ni in s and d 
band representation with hybridisation taken into 
considers cioij is given by. 


H = S ^ ( r ) P 

n ^ ns 


2 V® ( r -r ) T® t ^ r )P 
m^n ^ n mn n ^^n"^ ’ 

m^n 


nd 


2 


(r -r ) -!-2 y( r ) ( -t- ) 
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o ^ 

where P^^^ are projection operators and 

s d ds 

» '^n transfer operators defined by the follovdng 

expressions 


Pns ='ns><nsl j 


[ nd > <ndj 1 ^ T, 


with n,m = site indexj s,d = band 
indices. 


= ! ms > < nsl 

= 1 nd> < md 

= I ns > < nd 
= lnd> < ns 


rnn 

fd 

nm 

_sd 

^n 

,ds 


Here, site energy and hopping integral term (V) of 

the system can be considered as matrices in the following 
ways 


Si =sTi) 


/ 

/ 

( 


s / \ sd/ •* s \ 

( 3:,- ) - Y ( r. ) \ 


! sd/ -♦ 1 
\”Y ( ) 

\ 


i 

„d/ ■* s / 

\(rp ; 


(4.85) 



and 


■ / 




[v.j]= vcJ-.Jj) 


^ irCiS/ 

\ V ( Ti-rj ; 


v^V?,) \ 


V^Ci^.r.) I 


/ (4,86) 

Note that, if one is interested in writing the d band part 
of the Hamiltonian in equations ,(3, 7, 3, 8) for its two d- 
components ( t 2 g» ©g) along with its s-band part, then 
would be a matrix 3x3 in s, e^t ^2g ^representations. 

Here probability of occupation of site energy e^ 
is, given by,. 


Pr(ej_) = ('•-c)6(ej^ - e^) 


( 


where, 




S>cl , S^ci \ 

= e^’ n. T- (l-n.) 


'B 


and 


'"’‘’{Ti-r.) = n.nj V|>'1 (1-np(1-nj) 

+ ''ab'* ^ 


V 


s , d _j_ .^^s , d| 
BB 


+ V^*^(n^+n^) + V|»^ n^n 


where 


V 


(4 .$9) 


= 2 Vis'* 


and n- or n . represents random occupation variables having 

J 

values 0,1 for binary alloys. 



For the calculate onal part we have considered 
cluster effects only from randomness in the d-band and 
treated s-band and s-d hybridisation in CPA (since the 
off-diagonal disorder corresponding to s-band is very 
small) , 

We start from the Kubo-Greenwood form of the 
conductivity as before, 


= / dE (- ||) /fE) (4.91 ) 

-oo 

where,- f(E) = Fermi dirac di s tiibuti on, oc , P being 
directional dependences relative to the field direction* 


Now 


ap 


2“^E) 


Total = t scf'*‘°ds^ (4.92) 

conductivity 


= 2m Tr li“ 6 (EI-H) (El . H) ) 


= C|) Re Tr G (e") / G(E") - rgCE"^) /g(E'^)} 

es* CS>€%> 


{ S“P(E^, E“) - S (e”^, e"^) } 


(4,93) 


aS / % cc 

S is a ( 2 x 2 ) matrix in the basis s and d through i 

and (4 matrices. 


In general if one is interested in 
component of density of states or conductivity then S 
would be 3x3 matrix in the basis s, Sg and ^ 2 g * ^ 

particular matrix representation of S is given by 

= I Re /(d^Ic /8 7^) ^(k) ( 4.,94) 



Following the earlier equations (4,60 & 4,61 ) the Ward identity 
and Bethe-Salpeter equation corresponding to a multiband 
system has the following forms 

, { y. . z„sk) = -f S 

lj.y 


5 z (z.^z^fk) = f , , (.z.,z -Xk’) 

^ ^ Sn'^ ^ T 4 


^ a^' ^ > 22 J k* ) (4 ,95 ) 


ind 


( 1 ) 
V V 


b. ' (z^ , Z2I k) 


^ ( s ; - , ^ ( sXs ) , 

= S j (k) A ( k) +2 2 A (z. ,z^jk) 

a,p ar] \ir\^av ^ ria ^ p.p'aa' ' 


^hcT»«6 (z^,Z2|k ,k’) (z^,Z2|k) 

(4.96) 

All the indices a r r] , \i \ v ' , 5 , d , v , etc. refer to the 
appropriate bands of the systemi 

Definition of A (2^,, Z2|k) is given by 


rPd ) 


Q(’i > >=); JgC ic ) S ( Z 2 ! k ) 

. we may convert ■ super-matrices to matrices 


(4.97) 


As before 
and matrices to vectors as follows 


then 


«Sp = SSds' ®°dd) 

( 6 ) ^ 

(z^,Z2;k) = |^(d^k’ )/8'f^) Ap^( , Z2I k, k' ) 

6 g (^. 99 ) 



Now, equations for vertex furrction L have the following 
forms 


( 1 ) 

z -k) (z ^ 4) + 

( s ) -♦ d^k’ t -* 

^ qr "^pq^' ^1 ’ ^ ^'qr ^ ' ^2* ^ ( z^ , z^^k' ) 

P O) 

where L° = First term in the equation for \ (4,100) 


Here, the notations p, q, r refer to the super matrix 
indices which can take ss, sd, ds, and dd notations at a 
time as it is clear from representation of A in the 


following equation; 


A ■ 

A 

ss, ss 

ss, sd 

A , 

A 

sd, ss 

sd, sd 

A , 

A 

ds, ss 

ds, sd 

A , , 

A 

dd, ss 

dd, sd 


ss, ds 

A o , 

ss, dd 

^'^sd, ds 

^sd, dd 

^ds , ds 

A 

ds, dd 

^dd, ds 

A ' 

" dd, dd ! 


(4,101) 


Note that Lp, have also 
likeds t'Sgp in s and d Land 
The various Green Functions 

/ 2CP 

(k|E) = I Gg^ 


i 2CP 



four-vector representation 
indices in equation (4.98 and 4^99), 


can be evaluated as follows; 


2CP \ 

ii 

G 




I 


s( k) , 


sd 

2CP 


: = ■■li 


^sd 


\ 


-1 


I “^ds 


z-J)^-S,ds( k) / 


G 


/ 


dd I 


(4.102) 



where 


(k,E) 


„ d „ d , 7 V, 

z — S ^ S ^ s( k) 


d 


[{z -S^ s(k)}{.z-2, -EihS(10}-S,h 3 


P 2CP f V . p 'j 
Gdd 


z - S® - S® s(k) 


1d' 


^sd 


( 4.103) 


[!z - 4 - SijsC'WJiz-i;^ - 


(4.104) 


GlfihE) = —5 — ^ 

[{2-S® - 2^5s(k)}{z-S^ -s^^s(k)} ] 


sd 
(4.1,05) 


Now, 2CPA Green Functions corresponding to s,d band 
indices can be written as the following integral 
representations? 

-I- 00 


+ 00 


G^^(E) = G^^ (T1,E) Pg(Tl) dlj G5^^(E) = / G^-^(n,E) 

and d 


.2CP 


.2CP/ 


+ CO 


( n ) d t ) 


.2CP 


■00 


(4.106) 


where, 


Gg3^(Tl?E) 


Z -2^-Sid(n-j5/p3 


( Z-Tl) { Z - S 




^1 




z -S o “ 2^g(T)- SQ)/2^d 


(z -r|) iz - S' 


2is(n- s° ) 


'Id 


.j (4.107) 


(4.108) 


The derivation of the above results have been shown in Chapter 3 



Gg^^CniE) 


sd 




(4,109) 


We note, as before, (Mukherjee et al 1985, Rubio 1964) that the 
study of scattering diagrams for the kernel in equation 
(4,95) indie a tos that it is separable. So, we have. 


Apq ( z,, , 22? k, k' ) = ( z,, , Tc) a^ (z2,k») 

,+ 

Z2 


( 4 . 110 ) 


Now if 2-[ = E 


(^) ds (E,k) / 

ap= aO) (E,k) = / 




d\t %(E»k’) dgp(E,k’) 
S / dE dE"^ 


q 8 


( 4 . 111 ) 


And if z,| 



Note that from 

M 

. represents 
If 1 

band indices 
basis for the c 


Z2 = E”, Sp - a^^^ (E,k) 

Im 3^^ ( E, k) 

' '" (" s ' J" " g" " — 

g (E,k«)Im gJE,k')>'' 
q 8 ^ P ^ 

equation (4,95) to equation (4,112) notation 
summation and I'p,!; represents self-energies w.r.t, 

, etc. Equations (4,100) to (4,112) form the ; 

clculution of conductivity through the following 


steps? 

(i) Calculating the vertex correction from the Ward 
Identity in Equation (4.95) using CCPA self-energies? 

(ii) Iterative solution for vertex function L using 
vertex corractic«a frcwi (i)j 



(iii) Then using vertex function L to calculate 
the correlation function S and hence o(E), 

Results and Discussion on Static Conductivity Results 
xn Cu^_^Ni^; 

It has been pointed out by various authors (Niizeki 
1977a,b,c| Roth and Singh 1982) the vertex correction is a 
non-vanishing quantity in disordered' systems when we 
do conductivity calculations* We have pursued the calcu- 
lation in general in the light of multiple-band version 
of cluster CPA, In the previous section, we stated that 
CPA density of states for Cu. Ni (with hybridisation 
effects) is featureless. That is also true for conductivity 
in CPA as shown in Figure (4^7a) . As we go over to 
cluster CPA, the vertex correction is non-vanishing and 
negative for energies as it is related to backsca ttering. 

Also, the overlap terms between Cu and Mi are not sufficiently 
different to give an appreciable current correction, v\tiich 
is an order less than the vertex correction ^ 

?ife have calculated the 2CPA conductivity with and 
without correction for the concentrations c = 10%, 50%^-f 4-7 b»C) 
For calculational purposes we have solved the d-band part 
of the self-energy, i,e,2^, S^^in 2CPA while s-band and 
s—d hybridisation has been treated in single site CPA, 

As stated in Chapter 3, s-band is featureless in the pure 
system (Cu) and the off-diagcnal disorder corresponding to 




Fig. 4*7 a Static « 
CPA as 






Energy in Rydbergs 

c Static conductivity cr'CE) in the impurity band regime 
for c s 50 V« Nl In CU|_^ Ni^ olloy 
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s-band in the constituents Cu, Ni is less. 

The effect of vertex correction and statistical 
clustering both of them play important role in producing 
sharp features in conductivity. One can see from the 
Figure (4 ,7c, d) that almost all the features of density of 
states in 2CPA is reflected in the conductivity 
except in a small region where the vertex correction seems 
to be large enough to produce a dip at c = 50%, At c •* 1 0% 
(See Figure 4,7b ) the vertex correction is though small 

but the cluster effect which comes from the self-energies 
via Green Function is prominent to produce most of the 
structures in the conductivity. Hence the cluster effect 

in the conductivity formulation comes through self-energies 

di 

vertex correction should take care 

of them atleast doing 2CPA, We are able to yield analytic 
results for each of these concentrations throughout the 
band. 

Conclusions ? 

Now the calculation of electrical conductivity 
in CCPA yields the following interesting results: 

(i) Most of the fine structures in 2CPA density 
of states for c = 10% and c = 50% is reflected in 
conductivity without vertex-correction as shavn in Figure 



139 


s-band in the constituents Cu, Ni is less* 

The effect of vertex correction and statistical 

clustering both of them play important role in producing 
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Figure (4 ,7 c,d) that almost all the features of density of 

states in 2CPA is reflected in the conductivity 

except in a small region where the vertex correction seems 

to be large enough to produce a dip at c = 5C^, At c 1 0% 

(See Figure 4,7b ) the vertex correction is though small 

but the cluster effect which comes from the self-energies 

via Green Function is prominent to produce most of the 

structures in the conductivity. Hence the cluster effect 

in the conductivity formulation comes through self-energies 
d 

vertex correction ^^^^. 2 * 0 ne should take care 
of them atleast doing 2CPA, We are able to yield analytic 
results for each of these concentrations throughout the 
band. 

Conclusions ; 

Nov/ the calculation of electrical conductivity 
in CCPA yields the follov/ing interesting results; 

(i) Most of the fine structures in 2CPA density 
of states for c = 1 0% and c = 50/^ is reflected in 
conductivity without vertex-correction as shown in Figure 



All the previous cluster CPA calculations were done in 
alloys having featureless s-bands* Now we are in a stage to 
compute electronic density of states v^^ithin CCPA for(Cu^_^Ni^) 
realistic band like d-band which has considerable structure 
in the pure case. The density of states in the impurity 
band shows fine structures resembling the effects of off- 
diagcnal disorder and statistical clustering. 

Regarding the fine structure in electronic density 
of states of random alloys Faulkner (1982) made the following 
remarks, "Experiments to measure the density of states n(E) 
in an alloy all have a finite resolving power, and will 
therefore yield a smoothed average of the actual n(E), 

Since the CPA yields an excellent smoothed average of the 
density of states, it is expected to agree well with 
experiment and this does not mean that the resonance structure 
is unimportant". According to Faulkner ( 1 982) , the electronic 
driving force that would cause an order- disorder transfor- 
mation in a crystal has its origin in the resonance structure 
in the density of states, 

(iii) The earlier conductivity calculations 
( Brouers and Vedyayeb 1972, Niizeki 1 977a, b,c, and Niizeki 
& Hoshino 1977 ) are not of much realistic value because 
of the use of very crude model density of states in Brouers 
et.al.'s case or in Niizeki' s work with the choice of bond 
CPA which does not represent any realistic system results. 

Our calculation is realistic because we have used firstly the 



CHAPTER 5 


CONCLUDING REI^RKS 


In this work we have carried out the application of 
Augmented Space Formalism (ASF) to electronic density of 
states and static conductivity for binary alloys,. 

Although the feasibility of ASF has been previously 
proved for the calculation of electronic density of states 
in various model systems (See Chapters 2 and 3), we think 
that its feasibility for producing analytic results for 
the same for realistic models of metallic alloys like 
Cu^__^Ni^ provides a new way of understanding the underlying-, 
physics of cluster induced states in metallic alloys,. 

Here, we have performed the calculation of electronic 
density of states for a s— d band modeJ. Hamiltonian with the 
realistic choice of the various tight binding band parameters 
and using a 10 level of Recursion coefficients for the 
d-band for a finite size of cluster of f,,c,c. lattice,., 

The spiky nature of density of states is sraoothened by 

using the super-smooth program, . provided by Cambridge Recursion 

Library ( Solid State Physics, 35, 1980, p,78 , C,M,M. 

Nex 1978, Haydock and Nex 1984). This is facilitated 
by the significant numerical development made recently 
in. Recursion Methods. (Haydock and Nex 1985) through devising 
general terminator scheme for the Green Functions suitably for 



realistic d-band density of states for the host (Cu), 
secondly we have investigated the effect of vertex correction 
(which comes through two-site CPA self-energies) in the 
conductivity going beyond the traditional CPA theories 
(i,e, like single-site CPA theory or the CPA by Blackman 
et al, 1971) to consider off-diagonal disorder within ASF, 
Conductivity is sensitive to the effects of randomness, 
in current as well as to the vertex correction for strong 
of f— diag Qua 1 disorder which has been proved for the model 
simple cubic lattice of a disordered alloy. In the Cu. Ni 

! •*' C C 

disorder alloy vertex correction produces more significant 
effects in the conductivity compared to the effect of random- 
ness in current. 

We hope that the above investigations would thxow 
some lii^t in understanding the physics of electronic 
structure and transport properties of realistic alloys. 



Numerical Problfe m in the Calculation of DensitY_.^f 
a.i^---^-;^^~.~^-onduc.j.j_vj_-t^y within ASF t 

We think that one should study electronic density of 
states and static condictivity as an energy functional for 
larger cluster sizes within ASF aa the future plan. We have 
studied cluster effects at its minimum level. The problem 
becomes highly cumbersome if one increases the size of the 
cluster as this takes very long computation time to yield 
self-consistent solution for self— energies corresponding to 
a set of coupled nonlinear equations. In our calculation^ 
the average elapsed time for generating density of states 
over the entire band for a single concentration is approxi- 
mately one hour which is reasonable. The calculations may 
be considerably reduced by using real space symmetries of 
the cluster. This has been attempted so far as we know 
only in recursion .The future extension must include 
developments along this line if large clustereffects are to 
be made feasible. 

The main difficulty of solving vertex correction 
and finally the conductivity lies in the evaluation of 
multiple Brilouine Zone ( BZ) integration which appear in 
equations (4 .60,4 .61v4-.7l)For a system like Cu^ J'Ji^ where 
Feirai-surface can be considered approximately as a sphere 
the BZ integrals can be simplified considerably. But for 
many systems this may not be true. Note that, since here 
the integrand for these BZ integrals does not have cubic 



symmetry, one can not apply special direction techniques 
(R, Prasad and A. Bansil, 1980) or even the ’Ray* integration 
(Chen, 1977) must be generalised. This poses a severe 
problem of solving vertex correction equations for general 
disorder alloys until such theoretical developments are 
carried out, * 

Although we have given here the general multiband 
formulation of cluster CPA density of states and electrical 
conductivity yet for simplicity we have considered s-band 
and hybridisation effects in the spirit of CPA while d-band 
has been treated fully (diagcnal and off-diagcnal disorder 
both) in 2CPA, 

Finally we conclude with the following statements 
about the achievements that have been made here in the 
theories of disordered alloys; 

(i) Cluster extension of coherent potential 
approximation (CCPA) within augmented space formulation 
(ASF) for disordered alloys is (a) Herglotz at all energies, 
concentrations and degrees of disorder, and (b) it is 
feasible for realistic calculations, 

(ii) The use of supersmooth library routines 
(Haydock and Hex 1984) have eliminated the spurious nature 
of the electronic density of states (which comes from the 
tiuncation of the Green Function at a finite level) and 
produces most of the structure in pure density of states. 
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both single and multibands. 

The progress of electronic transport theories in 
disordered alloys seems to be less compared to the theories 
in electronic density of states. Until now nobody made an 
attempt to discuss electronic transport (static conductivity, 
Hall effect) going beyond single-site etc, or single-band 
CPA theory to study the general features which may appear 
due to statistical clustering, short- ranged order and off- 
diagonal disorder etc, except in the other extreme regime of 
localised state using the "maximally crossed" diagrams 
(AALR 1979 ), The nonvanishing nature of the vertex 
correction (due to back-scattering effects) has been recognised 
by Roth (1974), Roth and Singh (19B'2), Niizeki (1977a, b,c,v 
Niizeki and Hoshino 1977) in electrical conductivity of 
alloys. Following this we have established its importance 
for a model system (simple cubic lattice) with strong 
off-diagonal disorder and finally for a realistic alloy 

CUi Ni in the frame work of cluster CPA theory within 

I •— c c 

ASF., 

Explicit representation of vertex correction as a 
function of self-energies and cluster CPA Green Functions 
has bem derived in Chapter 4, The calculation of these 
quantities in the multiband form within cluster CPA for 
otiher real systems (like Cu-|__^Pd^» ^'^l-c^c^ where off 
diagonal disorder is strong may lead to yield interesting 
results. This is a future program of research,^ 
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(ii) With the introduction of vertex correction 
the conductivity is diminished having dips very prominen’t 
in the impurity band regime showing deviation from 
structure of density of states as reflected in case (i)i 
So this says that vertex correction which is always negative 
for energies considered in our problem has very prominent 
effects on the conductivity. In realistic systems like 

‘^'^1-c^'^^c* ‘^'^1_c^*^c’ ^’^1-c^^c off-diagonal disorder 

and pair scattering play significant role in self- energies-, 
the vertex correction is also significant there. In the 
Cui system overlap terms between Cu and Ni are not 

sufficiently different to give an appreciable current 
correction which is an order less than the vertex correction. 
In systems Cu-|_^Au^, Cu^_^Rh^, Cu^_^Pd^ one should take care 
of both vertex correction and current correction atleast in 
the spirit of 2CPA. 
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